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Abstract 

This paper analyses the GW method for finite electronic systems. In a first step, we pro¬ 
vide a mathematical framework for the usual one-body operators that appear naturally in 
many-body perturbation theory. We then discuss the GW equations which construct an ap¬ 
proximation of the one-body Green’s function, and give a rigorous mathematical formulation 
of these equations. Finally, we study the well-posedness of the GW° equations, proving the 
existence of a unique solution to these equations in a perturbative regime. 
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1 Introduction 


Computational quantum chemistry is nowadays a standard tool to numerically determine the 
properties of molecules. The Density Functional Theory (DFT) first developed by Hohenberg and 
Kohn [20] and by Kohn and Sham [25] , is a very powerful method to obtain ground state properties 
of molecular systems. However, it does not allow one to compute optical properties and electronic 
excited energies. In order to calculate such quantities, several approaches have been considered in 
the last decades [31]. Among them are the time-dependent DFT (TDDFT) [28, 29], wave-function 
methods [19] such as Coupled-Cluster, full-CI and Green’s function methods. In this article, we 
study the GW method, which is based on Hedin’s equations for the one-body Green’s function [18]. 
The formal derivation of the latter equations relies on many-body perturbation techniques. While 
the GW method has been proven very successful in practice to predict electronic-excited ener¬ 
gies, no rigorous mathematical framework has yet been developed to understand its mathematical 
properties. The aim of this work is to present such a framework. 

In non-relativistic first-principle molecular simulation, the electrons of a molecular system 
are described by an Wbody Hamiltonian operator H^, which is a bounded below self-adjoint 
operator on the fermionic space /\^ (see Equation (30) below). Whenever N ^ Z, where 

Z is the total nuclear charge of the molecular system, has an infinity of discrete eigenvalues 
E'^ < EIj < Ejj < • • • below the bottom of the essential spectrum, where E^ is its ground state 
energy. The quantities we would like to evaluate are the electronic-excitation energies 

E% — (gain of an electron), and E% — E%_^ (loss of an electron). 

These energy differences are not to be confused with the optical-excitation energies, which are en¬ 
ergy differences of the form E^ — E^^ between two states with the same number of electrons. More 
generally, it is interesting to compute the particle electronic-excitation set Sp := a [Hn+i — E^) 
and the hole electronic-excitation set S'h := cr — Hn-i)- As will be made clear in Section 3.2, 
these sets are closely linked to the one-body Green’s function: the time-Fourier transform of the 
Green’s function becomes singular on these sets. In order to study the electronic-excitation sets, 
we therefore study the one-body Green’s function. Also, the one-body Green’s function is a fun¬ 
damental object which contains a lot of useful information, and allows one to easily compute the 
ground state electronic density, the ground state one-body density matrix, and even the ground 
state energy thanks to the Galitskii-Migdal formula [15]. 

Galculating the one-body Green’s function is however a difficult task. In his pioneering work 
in 1965, Hedin proved that the Green’s function satisfies a set of (self-consistent) equations, now 
called the Hedin’s equations [18]. These equations link many operator-valued distributions, namely 
the reducible and irreducible polarizability operators, the dynamically screened interaction oper¬ 
ator, the self-energy operator, the vertex operator, and of course the one-body Green’s function. 
The state-of-the-art method to compute the one-body Green’s function consists in solving Hedin’s 
equations. 

Immediately, two difficulties arise. The first one is related to the lack of regularity of the Green’s 
function (we expect its time-Fourier transform G to be singular on the electronic-excitation sets). 
One way to get around this problem is to consider the analytical extension of G into the complex 
plane, which we denote by G. This is possible whenever the following classical stability condition 
holds true^: 


Stability assumption: It holds that 2E% < E%_pi E%_i 


^The question “Is the stability condition always true for Coulomb systems” is still an open problem [3, Part 
VII]. 
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The physical relevance of this inequality is discussed for instance in [11, Section 4.2]. It allows one 
to define the chemical potential /r, chosen such that 

~ ^N-l < M < ^^+1 ~ ^N- 

Instead of studying the Green’s function G(t) in the time domain, or its Fourier transform G(w) 
in the frequency domain, we rather study its analytical continuation G on the imaginary axis 
/i + iK. The function uj i—>■ G(/i + iw) enjoys very nice properties, both in terms of regularity and 
integrability, which makes it a privileged tool for numerical calculations. 

The second difficulty comes from the fact that Hedin’s equations cannot be exactly solved 
and, even more importantly, that the mathematical definition of some terms in these equations 
are unclear. It however opens the way to some approximate resolutions. The most widely used 
approximation nowadays is the so-called GW-approximation, also introduced by Hedin [18]. These 
equations are traditionally set on the time-axis, or on the energy-axis [37, 26]. However, as pre¬ 
viously mentioned, the various operators under consideration are singular on these axes, which 
makes the traditional GW equations cumbersome to implement numerically, and difficult to an¬ 
alyze mathematically. In order to manipulate better-behaved equations, it is more convenient to 
replace every operator-valued distribution involved in the GW equations by its analytic continua¬ 
tion on an appropriate imaginary axis, thanks to the “contour deformation” technique introduced 
in [36, 34]. The resulting GW equations, which give an approximation of the map w >->■ + 

turn out to give simulation results in very good agreement with experimental data [41, 42, 7, 8]. 

From the GW equations set on the imaginary axis, several further approximation may be per¬ 
formed. The GW equations are solved self-consistently, and the Green’s function is updated at 
each iteration until convergence. When only one iteration is performed, we obtain the one-shot 
GW approximation, also called the GqW® approximation of the Green’s function. For molecules, 
self-consistent GW approaches give results of similar quality as GqW®, sometimes almost identi¬ 
cal [42, 26], sometimes slightly worse [37], sometimes slightly better [7, 8]. When several iterations 
are performed, while keeping the screening operator W fixed, equal to a reference screening oper¬ 
ator VF°, we obtain the GW° approximation of the Green’s function [42, 46]. Since the update of 
the screening operator IF in a self-consistent GW scheme seems difficult to analyze mathemati¬ 
cally, we prefer to study in this article the equations resulting from the GW*^ approximation. 

The purpose of this article is threefold. First, we clarify the mathematical definitions and 
properties of the usual one-body operators involved in many-body perturbation theory. Then, we 
embed the GW° equations in a mathematical framework. Finally, we prove that, in a perturbative 
regime, the GW*^ equations admit a unique solution close to a reference Green’s function. 

From a physical viewpoint, the analysis we perform in this work is more relevant for atoms and 
molecules. Indeed, as discussed in [5, Section 4.1] for instance, fully self-consistent GW approaches 
are questionable for solid-state systems, for which quasiparticle methods are preferred [1, 2]. 

The paper is organized as follows. In Section 2, we provide the mathematical tools that will be 
used throughout the article. We recall the Titchmarsh’s theorem, and introduce the kernel-product 
of two operators, which can be seen as an infinite dimensional version of the Hadamard product 
for matrices. We also explain the underlying structure that makes the “contour deformation” 
possible. In Section 3, we recall the standard definitions of the usual one-body operators that 
appear in many-body perturbation theory. A consistent functional setting is given for each of 
these operators, and their basic properties are recalled and proved. Section 4 is concerned with 
the GW approximation. We explain why some of the GW equations are not well-understood 
mathematically, and prove that the GW° equations are well-posed in a perturbative regime. Most 
of the proofs are postponed until Section 6. 
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2 Setting the stage 


2.1 Some notation 


The GW method is based on time-dependent perturbation theory and therefore involves space- 
time operators. Following the common notation in physics, we denote by t the time coordinate, 
by r the space coordinates, and by x or rt the space-time coordinates. The functional spaces 
considered in this work are by default composed of complex-valued functions, unless we explicitly 
mention that the functions are real-valued. 

Most of the space-time operators appearing in the GW formalism are time-translation invariant. 
A time-translation invariant operator C can be characterized by the family of operators (C'(T))reR 
such that, formally, the kernel of C is of the form C(riti, r 2 t 2 ) = ^(ri, r2, where C(r, r', r) 

is the kernel of the operator C{t). For clarity, we will systematically use the letter r to denote a 
time variable which is in fact a time difference. 

Let 'H be a separable complex Hilbert space, whose associated scalar product is simply denoted 
by (•, •) and the associated norm || • ||. We denote by B{'H) the space of bounded linear operators 
on H, by S{'H) the space of bounded self-adjoint operators on H, by &p{H) (1 < p < oo) the 
Schatten class 


6pCH) = {a G B{'H) 


||A||e,(H) :=Tr(|Ani/^'<oo}, 


and by A* the adjoint of a linear operator A on T-L with dense domain. The real and imaginary 
parts of an operator A G B{'H) are defined as 


Re A = 


A + A* 
2 


Im A = 


A-A* 
2i 


Note that, when A is closed (which implies A** = A), the operators Re A and Im A are self-adjoint. 
For f,gG'H and given operators A, B on 'H, we will often use the notation 


{f\A\9)n ■■= if, A9)n, {f\AB\g)n := (/, ABg)n, 


even in cases when the operators are not self-adjoint. Operators are always understood to act on 
the function on the right in this notation. 

We will sometimes need to manipulate the adjoints of operators between two different Hilbert 
spaces Ha and Ht- The adjoint of a bounded operator A G B{Ha,Hb) is the bounded operator 
A* G B{Hb,Ha) defined by 

V(x, y) GHaX Hb, {A*y, = {y, Ax)^^ . 

Let A be a Banach space. We denote by o5^'(IR, E) the space of A-valued tempered-distributions 
on R, i.e. the set of continuous linear maps from the Schwartz’s functional space into E. 

Recall that, by definition, a family (T^)^>o of elements of E) converges in E) to some 

T G A^'(R, E) when rj goes to 0 if 

v<^g.^(R), ^ 0. 

Let / G L^(R, FI) be a time-dependent A-valued integrable function. The time-Fourier trans¬ 
form of / is defined, using the standard convention in physics, as 

Vw G R, /(w) := {Erf) (w) := [ dr. (1) 

Jr 


For the sake of clarity, we will sometimes denote by Rt or R^ the time-domain, by R^j the frequency- 
domain, by c5^'(Rr, £^) (resp. E)) the space of time-dependent (resp. frequency-dependent) 

A-valued distributions, etc. We will also denote with a hat the functions defined on the frequency 
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domain. Using this notation, can be extended to a bicontinuous isomorphism from E) 

into When / G we have 


Vr G 


Et^T) (r) 


i 


/(w)e-'“^ dw. 


The Dirac distribution at a G is denoted by Sa, and the Heaviside function on R by 0: 

0(r) = Ifor T > 0, 0('r) = 0 for r < 0, 0(0) = 1/2. (2) 

Recall that the time-Fourier transform of 0 is, in the tempered distributional sense, 


0(a;) = ttSo{uj) + ip.v. ( ^ ) , 


(3) 


where p.v. is the Cauchy principal value. We will also make use of the notation r'*' for a number 
strictly above t, but infinitesimally close to t, and of the convention 

0( r )(5 o ( r +) := So{t), 0(- T )(5 o ( r +) := 0. 

2.2 Hilbert transform of functions and distributions 

The Hilbert transform, which amounts to a convolution by 7r“^p.v.(i), plays a crucial role in 
the GW formalism. We first recall some well-known results on the standard Hilbert transform 
on LP(R^), and extend the results to the Sobolev spaces it®(R(^) for s G R. Usually, the name 
“Hilbert transform” is only used on functional spaces E C such that, for any function 

f G E, the limit 


/*p.v. ( - 


J-oo ^ ^ ri^0+ J^\ 


/K) 


[ui-ri,uj+ri\ W — W' 


dw' 


exists for almost all w G R^j. However, in the sequel, we will also use the name “Hilbert transform” 
in functional spaces where the above integral representation is not always valid (for instance when 
/ is not a locally integrable function). Note that we define the Hilbert transform on Fourier 
transforms of functions (*.e. on functions on the frequency domain) since this is the typical 
setting in the GW formalism. 


2.2.1 Hilbert transform in spaces 

We first begin with the following classical definition (see for instance [17, Section 4.1]). 

Definition 1 (Hilbert transform on .5^(Rpj)). The Hilbert transform of a function (p G is 

defined by 

-p-v. (4) 


or equivalently by 


b</ := (^^(-isg^-)) p, 


(5) 


where p.v. (i) is the Cauchy principal value of the function a; * the convolution product, Tt 

the Fourier transform defined in (1) and —isgn(-) the multiplication operator by the L°° function 
1 1 —>■ —isgn(t) (where sgn(t) = 0(t) — 0(—t) is the sign function). 

The Hilbert transform can be extended by continuity to a large class of tempered distributions. 
We refer to [17, 35] for a proof of the following theorem. 
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Theorem 2. For all f G 


,) with 1 < p < oo, the Hilbert transform 


S)f{uj) = p.v 


f 




UJ — UJ' 


Alo' 


is well-defined for almost all w G K. It holds io S with 


tan(7r/(2p)) if 1 < p Gi 2, 
cotan(7r/(2p)) if2^p<oo. 


Moreover, the Hilbert transform commutes with the translations and the positive dilations, and 
anticommutes with the reflexions. Finally, it is a unitary operator on 


2.2.2 Hilbert transform in Sobolev spaces 

Recall that for any s G R, the Sobolev space is the Hilbert space defined as 

H%R^) := {/G y'{R^)) I (1 + I • G , 

and endowed with the scalar product 

^ + CXD _ 


/.?)^^ =27r/ (1 + r /)(r) (Tj. ?)(r)dT, 


and that can be identified with the dual of H’^iR,^) when the space L'^{Rui) = H'^{Roj) 

is used as a pivoting space. One of the reasons to introduce these spaces is that the image of 
L°°(Rt) by the Fourier transform Hr is contained in the Sobolev spaces of indices strictly lower 
than —1/2. 

Lemma 3 (Fourier transform in L°°(]Rt-)). Let s > 1/2. Then Ft C H~^(R,^) and 

( f cIt” \ 1/2 

ll■^T|lB(L-.ff-o) =C'« = (^27r • (6) 

For completeness, we recall the proof of Lemma 3 in Section 6.1. 

Since the Hilbert transform in S^(R,^) amounts to a multiplication by the bounded function 
—isgn(-) in the time domain (see (5)), it can be directly extended to the Sobolev spaces iL®(IR^). 

Lemma 4. For any s G R, the Hilbert transform f) is a unitary operator on the Sobolev spaces H’^iRF) 
satisyfing = —io (and therefore = —Id/. 

Remark 5 (Hilbert transform of distributions). Extending the Hilbert transform to Sobolev spaces 
is straightforward using (5). Extensions of the Hilbert transform to other subspaces of ^'(Ri^), 
such as the spaces defined in [38, Section VI.8], can be obtained from ///. 


2.2.3 Hilbert transforms of operator-valued distributions 


We now need to properly define the Hilbert transform of operator-valued distributions on the 
frequency domain, as such objects naturally appear in the GW formalism. We first introduce, for 
s G R, the Banach space 


HyR^,Bin)) := {1 G y'{R^,Bin))) I (1 + I • \y/^Ffi^A G LyRr,B{n))], 

endowed with the norm 

2 


A 




V^(^£°°il+ry\\{Ffi^A){r) 


Bin) 


dr 


1/2 


The following definition makes sense in view of Lemma 4. 
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Definition 6 (Hilbert transforms of frequency-dependent operators). Let % he a Hilhert space, 
and eonsider s G M and A G S('H)). The Hilhert transform of A is the element of 

denoted hy Sj{A), and defined hy 


y{f,g)enxn, 




9 = 


»((/ 



(7) 


In particular, it is possible to define the Hilbert transform of the Fourier transform of a uni¬ 
formly bounded field of time-dependent operators, using the following result, which is a straight¬ 
forward extension of Lemma 3. 


Lemma 7. Let TL he a Hilhert space, and let s > 1/2. Then for all A G L°°{M.r,B{'H)), we have 
A€ H-‘‘{R^,B{H)), with 




= 2tt 


r \ 

■Jjl + r^)-^\\A{r)\\ 

Bin) ^ II^I|l“(r^.b(h)) 


where Cg is defined in (6). 


Let ^(R) be the set of Borel subsets of R, 6 G i^(R.) a Borelian set, and H a self-adjoint 
operator on a Hilbert space H. We denote by := Ib(-ff) the spectral projection on 6 of TL 
(here, is the characteristic function of the set h, and tb{H) G B{'H) is defined by the spectral 
theorem for self-adjoint operators; see for instance [33, Theorem VH.2]). 


Definition 8 (Principal value of the resolvent of a self-adjoint operator). Let H he a self-adjoint 
operator on a Hilhert space TL. We define the B{Ti)-valued distribution p.v. the frequency 

domain by 


y{f,g)enxn, 



where is the finite complex Borel measure on R^^ defined by 


V&G^(R.), gfjb) = {f\P,^\g). 

As any complex-valued bounded Borel measure on R^^ is an element of iL“®(R^) for any s > 1/2 
(this is a consequence of the continuous embedding H^{Rgj) ^ C'°(Rtj) fl L°“(R) for s > 1/2), it 
follows from Definitions 6 and 8 that 


p.v. 



= 7rjo(P^) 


in H-^{R^,B{n)), 


s > 1/2, 


which is the operator analog of the well-known formula 


p.v. 



nS}{So) in H ®(Ra;), s > 1/2, 


( 8 ) 


which is itself a simple reformulation of the equality 





-^sgn(.) inL°°(R,). 


2.3 Causal and anti-causal operators 

The GW formalism makes use of families of time-dependent operators (Tc(r))rGR and (Ta(T)),-gR 
of the form 

Tc(t) = 0(r)Ac(T) and Ta(T) = 0(-T)Aa(T), 

where 0 : R —>■ R is the Heaviside function (2), and Ac and Aa belong to L°°{R,B{Tl)) for a 
given Hilbert space TL. The family of operators (T'c(T))r6R is called a causal operator, as Tc{t) = 0 
for all T < 0. Likewise, the family of operators (Ta(T)),-gR is called an anti-eausal operator, as 
Ta(T) = 0 for all r > 0. We recall in this section the basic properties of causal and anti-causal 
operators. 
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2.3.1 Causal operators 


Causal functions have very nice properties, because their Fourier transforms have analytic exten¬ 
sions in the upper half-plane 

U := {2 S C I Im z > 0} . 

This comes from the fact that, if / G L^(K.t) + T°°(K.r) is such that /(t) = 0 for r < 0, the 
Laplace transform / of /, defined on U by^ 

^ r r+oa 

VzeU, /(z) := / /(r)e‘"-dr= / /(r)e-" dr, (9) 

Jb Jq 

is a natural analytic lifting onto U of the time-Fourier transform / of / defined on = 9U. 
Note that the Laplace transform can be extended to appropriate classes of tempered distributions, 
see [38, Chapter VIII]. 

Let us first recall the Titchmarsh’s theorem [44] (see for instance [30, Section 1.6]). 

Theorem 9 (Titchmarsh’s theorem in [44]). Let f G L^(IR.t-) and f G L^(Mt^) be its time- 
Fourier transform. The following assertions are equivalent: 

(i) f is causal (i.e. fir) = 0 for almost all t < 0); 

(a) there exists an analytic function F in the upper half-plane U satisfying 

sup ( / 11^(0;-I-iry)]^ dw ) < 00 

r />0 \J-00 J 

and such that, F{- -|- iTy) —>■ / strongly in L^(]Rij), as g ^ 

(Hi) Re/ and Im/ satisfy the first Plemelj formula 

Ref = -Pj (im/) in 4v^(Rc.j); (10) 

(iv) Re/ and Im / satisfy the second Plemelj formula 

lm.f=Sj(Refj in L^(K,^). (11) 

If these four assertions are satisfied, then the function F in (ii) is unique, and coincides with the 
Laplace transform f of f. 

We refer to [44] for a proof of this theorem. Formulae (lO)-(ll) are sometimes referred to as 
the Kramers-Kronig formulae or the dispersion relations in the physics literature. Titchmarsh’s 
theorem implies in particular that square integrable causal functions, which can be very easily 
characterized in the time domain (they vanish for negative times), can also be easily characterized 
in the frequency domain (the imaginary parts of their Fourier transforms are the Hilbert transforms 
of their real parts). 

We emphasize that the above version of Titchmarsh’s theorem is only valid in while the 
GW setting mostly involves L°° causal functions (see Section 3.2 for instance). Weaker versions 
of Titchmarsh’s theorem are available for wider classes of tempered distributions (see [30] and 
references therein), but the setting turns out to be sufficient for our purposes and has the 
advantage of allowing short, self-contained proofs of all statements. Note that the assertions are 
no longer equivalent. 

^The Laplace transform is usually defined as 

POO 

F{p) = / /(r)e-^"dr. 

Our definition, which is better adapted to the GW framework, simply amounts to rotating the axis, or, in other 
words, to setting 2 = ip. 
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Theorem 10 (Titchmarsh’s theorem in L°°(R)). Let g G be a causal function (i.e. 

g{T) — 0 for T < 0) and let g G for all s > 1/2 be its time-Fourier transform, and g be 

its Laplace transform defined on U. Then, 

(i) g is analytic on U; 

(a) the function g g{- + ig) is continuous from ( 0 ,+oo) to for all s G R, and is 

uniformly continuous from [0,+oo) to il“®(Ri^) for all s > 1/2. Moreover, g{- + iry) —>■ 5 
strongly in for all s > 1/2, as g ^ O’*"; 

(Hi) for all z 

5(2) = ^ (?,(•-■ (12) 

(iv) Re 5 and Img satisfy the Plemelj formulae: 

Re^ = —.5 (Img) and ling = Sj(Reg) in il“^(R(j). (13) 


The proof of Theorem 10, which is a simplified version of the proof of the more general result 
given by [43, Lemma 1] (see also [30, Section 1.7]), is given in Section 6.2. For simplicity, we 
stated (12) and (13) in for the value s = 1, but similar results hold for any value s > ll2. 

Let us now extend these results to operator-valued functions. We recall that a map A{z) from 
an open set ?7 C C to a Banach space E is said to be strongly analytic on 17 if 17 9 z 1 —>■ A{z) G E 
is C-differentiable on U, i.e. 6.A{z)/<iz G E for all z G U. 


Definition 11 (bounded causal operator). Let TL he a Hilbert space and G L°°(Mt, B(H)). We 
say that Tc is a bounded causal operator on % if Tc{t) = 0 for almost all t < 0. 

Lemma 3 and Theorem 10 can be straightforwardly extended to operator-valued maps (see 
Section 6.3 for the proof). 

Proposition 12. Let H be a Hilbert space and Tc G L°°{M.t,B{'H)) a bounded causal operator 
on TL. Then its time-Fourier transform Tc belongs to H~‘‘(M.c,B{'H)) for any s > 1/2, and its 
Laplace transform 

^ p r+oc 

Tc{z) := / Tc{t) e'^" dr = / Tc{t) e‘^" dr 
7r 7o 

is well defined on the upper-half plane U. Moreover, 

(i) Tc is a strongly analytic function from U to B(H); 

(a) the function g 1 —>■ Tc(- -|- ig) is continuous from (0,-|-oo) to i7®(Rtj,S('H)) for all s G R, 
and uniformly continuous from [0,-|-oo) to H~‘^(E.,^,B{'H)) for s > 1/2. Moreover, for any 
s > 1/2, Tc{- -\- ig) —>■ Tc strongly in i7“®(R^, S('H)) as g ^ 0 +; 


(Hi) for all z G\], it holds 


Tc{z) 





(iv) the operators ReTc and ImTc satisfy the Plemelj formulae: 

Refc = -fi) (imfc) and Imfc = (Ref^) in H-\Rc,B{n)). (14) 


Besides the general case covered by Proposition 12, the particular case of causal time-propagators 
is often encountered. Explicit formulae can be provided for the Laplace and Fourier transforms in 
this case, as made precise in the following result (see Section 6.4 for the proof). 
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Proposition 13 (Analytic extension of causal time propagators). Let H be a self-adjoint operator 
on a Hilbert space H and Ac{t) := —i0(T)e“‘'^^. The Laplace transform (Ac(z)) 2 gu coincides with 
the resolvent of H in U; 

Ac{z) = {z- H)~^. 

Moreover, Ac{- + irf) converge to Ac in H~^as rj ^ 0+, and 
Re Ac = p.v. ^ImAc = —ttP^ in 

Let us conclude this section with a useful result (see Section 6.5 for the proof). 

Lemma 14. Consider a bounded eausal operator Tc S L°“(Rt,'B(’H)) sueh that Supp(ImTc) C 
[wo,oo) for some uiq € R. Then ImTc ^ 0 on if and only ifKeTc ^ 0 on (—oo,wo]. 


2.3.2 Anti-causal operators 

Definition 15 (bounded anti-causal operator). Let H be a Hilbert spaee and T^ S L°°{M.r,B{'H)). 
We say that Ta is a bounded anti-causal operator if Tg,{T) = 0 for almost all r > 0. 

All the results for causal operators stated in the previous section can be straightforwardly 
transposed to anti-causal operators, by remarking that if (Ta(T)),-gR is an anti-causal operator, 
then (Ta(—r))tgR is a causal operator. We will use in particular the following results, which are 
the counterparts of Proposition 12, Proposition 13 and Lemma 14. 

Proposition 16. LetH. be a Hilbert space and Tg, S L°“(]R.,-, S('H)) a bounded anti-causal operator 
on TL. Then its time-Fourier transform Ta belongs to H~^{]S.cj,B{'H)) for any s > 1/2, and its 
Laplaee transform Ta is well defined on the lower half-plane 

L = {2 G C I Im (z) < 0} . 


Moreover, 

(i) Tc is a strongly analytic function from L to B{TL); 

(ii) the funetion rj 1 —>■ Ta(- — irf) is continuous from ( 0 ,-l-oo) to H‘^(/E.,^,B{TL)) for all s G R, 
and uniformly continuous from [0,-|-oo) to H~‘^{M.c,B{'H)) for s > 1/2. Moreover, for any 
s > 1/2, Ta(' — ir]) — >■ Ta strongly in iL“®(Raj,yB('H)) as p ^ O'*'; 


(Hi) for a/Z 2 G L, it holds 

(iv) the operators ReTa and ImTa satisfy the Plemelj formulae: 


T.iz) = -^(Tc,f-z)-^) 

2 i 7 r \ /H-^,m 


Rei; =15 (imTk) and Imi; =-ii (Rei;) in H-\^c.,B{'H)). (15) 


Note that the signs in the Plemelj formulae are different for causal and anti-causal operators 
(compare (14) and (15)). Also, the Laplace transform is defined in the lower-half plane L for 
anti-causal operators, while it is defined in the upper-half plane U for causal operators. The 
counterpart of Proposition 13 is the following proposition. 

Proposition 17 (Analytic extension of anti-causal time propagators). Let H be a self-adjoint 
operator on a Hilbert space TL and Ac{t) := i0(—T)e''^'^. The Laplace transform (Aa( 2 :)) 2 gL Is 






Moreover, A^f — ip) converge to Aa in H ^(Rcj,B(/H)) as p ^ 0+, and 


ReAa = p.v. 


j—and ImAa = 7 rT ^ in H ^{M.cj,B{'H)). 
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Finally, a result similar to Lemma 14 can also be stated. 

Lemma 18. Consider a bounded anti-causal operatorTs, S such t/iat Supp(ImTa) C 

(— 00 , Wo] for some wq S Then, ImTa ^ 0 if and only ifKeTi,{uj) ^ 0 on [wo,+oo). 

2.4 Operators defined by kernel products 

Two of the fundamental equations in the GW method (see Sections 4.2.1 and 4.3.1) are of the 
form 

C(xi,X2) = i.4(xi,X2)-B(x2,xi), (16) 

where ^(x, x') and S(x, x') are the kernels of space-time operators invariant by time translations. 
As the product of the kernels of two operators is not, in general, the kernel of a well-defined 
operator, we have to clarify the meaning of (16). We first treat the case of time-independent 
operators in Section 2.4.1, and consider time-dependent operators and their Laplace transforms 
in a second step (see Section 2.4.3). 

2.4.1 Definition of the kernel product 

We first consider the special case when the operators in (16) are time-independent. Our aim is to 
give a meaning to equalities such as 


C(ri,r2) := A(ri,r2)B(r2,ri), (17) 

where A(r,r') and are the kernels of two integral operators A and B on L^(IR^). For this 

purpose, we replace (17) by the formally equivalent definition 

y{f,g) & X {f\c\g) := [ [ 7(ri)C'(ri, r 2 )g(r 2 ) dri dra 

= f [ A(ri,r2)5(r2)B(r2,ri)7(ri)dri dr2 

= Tri.(R3) {AgBj) , (18) 

where the last line involves the operators A and B themselves, and not their kernels (/ and g are 
there seen as multiplication operators by the functions / and g respectively). 

The formal equalities leading to (18) suggest to define the kernel product of two operators A and 
B (defined on dense subspaces of L^(M^)), as the operator on L^(R^) with domain L> C L^(]R^), 
denoted hy Aq B and characterized by 

V(/,g)GL2(M3)xD, {f\{AQB)\g)-.= TvL^R^){AgB'f). (19) 

In particular, the product A0i? is a well-defined bounded operator on L^(K^) as soon as AgBf is 
trace-class for all (/,<?) G L^(M^) x L^(]R^) and (/, g) !->■ Tr^2(R3)(Agi?/) is a continuous sesquilin- 
ear form on (R^) x (R^). ft follows from the above considerations that if A and B are operators 
with well-behaved (for instance smooth and compactly supported) kernels A(ri, r 2 ) and i3(ri, r 2 ), 
then A 0 i? is a bounded integral operator with kernel (A 0 B){ri,r2) = A(ri, r2)B(r2, ri). 

Remark 19. It is also possible to rely on the formal equality 

y{f,g) G L\R^) X l2(R3), {f\c\g) = Tri2(*3) (jAgB) , 

and define another kernel product 0 by 

V(/, 5 ) G L2(r3) X D, (/ |A0B| g) := Tv^r^) (jAgB) . 
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It may hold that AqB is a well-defined bounded operator, while AqB is an unbounded operator.^ 
In the sequel, we will mostly state the results for the 0 kernel product. 

Remark 20. The product AqB can be seen as an infinite-dimensional extension of the Hadamard 
product A o defined for two matrices A € ([ynxn g ^ £^nxm 

VI < i ^ m, (A o B^) = A^ (B"^) = A^ Bj*. 

Let us specify possible sufficient conditions for the operator A 0 _B to be well-defined. The 
typical situation we will encounter in the GW setting (see Sections 4.2.1 and 4.3.1) is the case 
when A € B(L^(R^)), while B is an operator on L^(K^) satisfying 

Vf,g€L^(R^), Ty{\gB7\) ^CsWfWML^- (20) 

In this case, the operator A Q B defined in (19) is a well-defined bounded linear operator on 
L2(R3), and 

m 0 B||s(L2(R3)) ^ ^311^110(^2(83)). 

The operators B arising in the GW formalism are usually of the form B = BIB 2 B 1 , where Bi 
is an operator from L^(R3) to some Hilbert space TL, and B 2 G In fact, assume that the 

operator Bi is such that Bif G ©2(L^(K.^),’H) for any / S L'^fR^), with 


\\Bif\\e2(L^(R-^).H) ^ K\\f\\L2, ( 21 ) 

for a constant K G independent of /. In the left-hand side of (21), / denotes the multiplication 
operator by the function /. In this case, (20) holds with 

Cb = K'^\\B2\\b(h)- 


Let us conclude by giving a simple example when (21) is satisfied, in the situation when 

n = l2(r3). 

Lemma 21. Let Bi be a linear operator with integral kernel Bi (r,r') e LL(M^ X R3)^ such 
that r H> IIHi(r, •)II 300 G L^(R^). Then Bi G B(L^(R^), L^(R^)), so that Bi defines an operator 
on L^(R3) yjith domain L^(R^) fl L^(R3). Moreover, for any f G L'^(R^), the operator Bif is 
Hilbert-Schmidt on L^(R3)^ with 


I|-Bi/||s2(L2(R3)) < 



The proof of this result can be read in Section 6.6. In the GW setting, a technical result similar 
to Lemma 21 is provided by Lemma 77. 


2.4.2 Properties of the kernel product 

Lemma 22. Consider two bounded operators A,B G B{L'^{R^)) such that A,B^0 and (20) 
holds. Then, Aq B is a bounded, positive operator on L^(R3^_ 

3 As an example of such a situation, take (f> S L3(r3) n L°°(R3), tp G L^(R3) \ and set A = IbKbl 

and B = \(f>){(j>\. Then, for all f,g € L^(R3), the operator AgBf = \'f){4’\g\4‘){SI\ S' well-defined rank-1 bounded 
operator since (pf S L^(R3), hence is trace class. Moreover, 

Tti2(K3) (AgBj) < (||blli~llblL 2 ||blL 2 ) ||/IL2||5iL2, 

so that A O S is a well-defined bounded operator on L^(R^). On the other hand, it formally holds fAgB = 
\fip){(l)\g\(l>){(l)\. If / is such that f'lp ^ then this operator is not bounded. 

We are grateful to Yanqi Qiu for pointing out this counter-example to our attention. 
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The proof of this result is very simple: it relies on the observation that, for any / S 
(/ 1^ QB\f) = Tri.(R3) (AfBj) = Tri.(R3) (^A^/^fBjA^/^) ^ 0, 

since fBf is a positive, trace class operator and A}!"^ > 0 is a bounded operator. 

Lemma 23. Consider two bounded operators A,Bg S(L^(K.^)) such that (20) holds. Then, AqB 
is a bounded operator with adjoint {A 0 B)* = A* 0 B*. 

The proof of this result is also elementary: for any f,gG 

(/| {AQB)g) = Tri2(R3) (AgBf) = Tri2(R3) (^{AgBf)*^ = Tri2(R3) {fB*gA*) 

= Tri2(R3) {A*fB*g) = {g,{A*QB*)f) = {{A* 0 B*)f, g). 

In particular, Aq B is self-adjoint whenever A and B are self-adjoint. 

2.4.3 Laplace transforms of kernel products 

We finally combine the results on causal operators with those on the kernel product 0 defined in 
Section 2.4.1 in order to give a meaning to (16). Note first that the space-time operator with kernel 
C(x, x') is also time-translation invariant and that the family of operators (24(T)),-gR, {B{T))r^m. 
and such that, formally, ^(xi,X2) = 24 (ri,r 2 ,ti - ^2), <B(xi,X2) = i?(ri,r2,ti - <2), 

and C(xi,X2) = C'(ri,r2,ti — t2), are related by 

Cir) =iAiT)QBi-T). ( 22 ) 

We assume here that A and B are such that (22) is well-defined. When all the operator-valued 
functions have sufficient regularity in time, their Fourier transforms decay sufficiently fast at 
infinity and it is possible to Fourier transform (22). This is however not the typical case in the 
GW setting since we work with causal and anti-causal operators, whose Fourier transforms are in 
iJ“®(IR^) for some s > 1/2. 

We therefore rather consider Laplace transforms. More precisely, for two helds of uniformly 
bounded operators (24(T))reR and {B{T))re»., and provided C'(r) := iA(r)0i?(—r) is well defined, 
we can decompose A, B and C as the sums of their causal and anti-causal parts as 

A(t) = A'^{t) A~{t) with A'^ij) := 0(T)A(r) and A~(r) := 0(—t)A(t), 

and similarly for B and C. Then, 


ij) = lA'^{ t) Q B (— t ) and C (r) = iA (T)QB'^(—r). (23) 


We next consider w > 0 and 0 < 77 < w. From the equality 


C+(t) e-“^ = i 


24 +( 


r e 




0 [B-(-T)e-^^] , 


we deduce, by Fourier transform, that 

__ i r+°° ~ __ 

C~^(i^ + iuj) = — / A+(ly — cj'+ i(uj — g)) Q B~ (—uj'— irj) dca'. (24) 

J —00 

The convolution on the right-hand side is well defined in view of Propositions 12 and 16. It 
however becomes ill-defined as 07,77 —>■ 0. In the case when the causal and anti-causal operators 
24+ and B~ under consideration are time-propagators, it is possible to remove this singularity by 
rewriting the convolution on appropriately shifted imaginary axes. 
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Theorem 24. Consider three Hilbert spaces and assume that 


A+{t) = 
B+{t) = 


A-{t) = 

B-(r) =i0(-T)S*e'^'®^Bi, 


where Ai G Bi G and ^2,-82 are possibly unbounded, self-adjoint operators 

on Ha and Hb respectively, for which there exist real numbers a, b such that A 2 ^ a and B 2 ^ b. 
We assume in addition that, for any f € H, Bif G & 2 (H,Hb) with ||i?i/||©2('H.'H6) ^ for 

a constant K G R’*' independent of f. Then, the operators C, C'^ and C~ in (22)-(23) are well- 
defined, the Laplace transforms of C'^ and C~ admit analytical continuations on 1[JULU(—00, a+6) 
and U U L U (—(a + b), 00 ) respectively, and it holds for any v < a-\-b and v' G (— &, a — v), 

Vw G R, C'^{v-\-iuj) = — — / A+{y-\-v'-\-\{ijj-\-uj')) (■) B~ {v'iu') Au', (25) 

J —00 


while, for any v > —(a + 5) and v' G (—a — v, b), 

_ 1 r+°° _ _ 

Vw G R, C~{v-\-iuj) =—— / A~ [v-\-v'-\-i{uj-\-ui')) (■) B+{v'-\-\u)')Auj'. (26) 

J — C50 

Finally, the following equality holds provided b > 0 and a + 6 > 0; for any v G (—(a + 5), a + 6) 
and v' G (—6, b), 

~ 1 /’+°° ~ ~ 

Vw G R, C'(^ + ia;) = —— / A{u-\-u'-\-i{ujuj')) Q B{v'-\-luj') Aui'. (27) 

J — 00 


The proof of Theorem 24 can be read in Section 6.7. The choices of i/, v' ensure that the 
function uj' !->■ A+{y-\-v' -\-\{uj-\-uj')) is in L'p(R,^,B{H)) for anyp > 1, while, for any f,gGH, the 
function oj' !->• gB~{C + iw')/ is in LP(R,^,&i{H)) for any p > 1. Therefore, in view of (25), the 
functions !->■ C'+(i^ + iw) is in LP(R,j, 8(77)) for any p > 1. Similar results hold for w !->■ C~{y-\-iijj) 
and to ^ C(y -\- iw). 

Let us conclude this section by deducing interesting properties from the analytic continuation 
results given by Theorem 24 (see Section 6.8 for the proof). 


Corollary 25. Assume that the conditions of Theorem 2j hold. Then, 

Supp ^Im ( 7 +^ C [a + 6 , +00) , Im C+ ^ 0, 

Supp ^Im(7“^ C (—00, — (a + 6)], Im(7“ ^ 0, 

so that 

Supp flm Cj C R\ (— (a + 5), a + 6), Im C > 0. 

Moreover, 

C+ = Re(7+^0 on ( — 00, a+ 6), 
C“=ReC“^0 on ( — (a + 6), +00). 

In particular, (7 = ReC ^ 0 on ( — (a + 6), o + 5). 


(28) 


(29) 


2.5 Second quantization formalism 

We recall here the definitions of the main mathematical objects used in the second quantization 
formalism, which are used to define - at least formally - the kernels of the operators arising in the 
GW method. More details about the second quantization formalism can be found e.g. in [9]. 
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We consider a system of N electrons in Coulomb interaction subjected to a time-independent 
real-valued external potential Vext G -I-R). In order to study the response of the 

system when electrons are added or removed, we embed this iV-body problem in a more general 
framework where the number of electrons is not prescribed. We denote by 'Hi = L^(IR^,C) the 
one-electron state space (the spin variable is omitted for simplicity), by Hat = 'Hi the N- 
electron state space, and by F = ©(^“qHiv the Fock space, with the convention that Ho = C. The 
Hamiltonian of the H-particle system reads 

1 W AT ^ 

Hn = “ 2 X/ X/ (’’i) + X/ Ij.. _ j. I > 

i=l i=l 

and the corresponding Hamiltonian acting on the Fock space is denoted by H, so that Hm = EI|wjv 


For / G Hi, the creation and annihilation operators a){f ) and a{f) are the bounded operators 
on the Fock space F defined by 

VH G N, a'^(/)|-H„ G H(HAr,HAr+i), a(/)|'H„^i G H(HAr+i,Hw), 
and for all ^at G Hat, 

, N+l 

[a1'(/)$Ar](ri,... riv+i) := ^== (-l)^’^V(i’i)«'Ar(ri,..., r^-i, r^+i,.. .,rN+i), 

^ + i=i (31) 

[a(/)d>iv](i’i, ■ •-rw-i) := [ /(r)$iv(r, ri,..., tat-i) dr. 

The creation and annihilation operators satisfy o^{f) = a{f)* and the anticommutation relations 

V(/,5 ) G Hi X Hi, [a{f ),aig)]+= 0, [a1'(/), a^(g)]+= 0, [a{f ),a'f ig)]+= {f\g)lw, (32) 

where \A, B\+ = AB -\- BA is the anti-commutator of the operators A and H, and where Ip is the 
identity operator on F. In particular, 

a\f)a(f) + a(f)a\f) = II/IIhi If- 

The mappings Hi 9 / !->■ a^(f) G H(F) and Hi 9 / !->■ a(f) G H(F) are respectively linear and 
antilinear. 

In most physics articles and textbooks, the GW formalism is presented in terms of the quantum 
field operators in the position representation d>(r) and 'l'l(r). We recall that, formally, 

OO OO 

VrGK^, 'i’Hr) = '^(j),(r)a'^((/),), ^'(r) = ^ ^i(r)a((/)i), 

2=1 2=1 


where is any orthonormal basis of "Hi. 

[ ^'^(i’)./’(r)dr = a'''(/) 

Jr3 

In the second-quantization formalism, HI reads, 


Note that for any f G Hi, 
and f 5'(r)/(r)dr = a{f). 

Jr3 


HI = 





■ 'Cext (l*. 


^'(r) dr -I- 


1 

2 


vl(t(r)vl''l'(r')|r — H| ^'I'(H)'I'(r) dr dr'. 


Finally, we introduce the Heisenberg representation of the annihilation and creation field op¬ 
erators 'I'H(rt) and 'l'Jj(rt), formally defined by 

5-tj(rt) =e™«'1'(i.)e-™ and «'H(rt) = e™«'(r)e-™. 

Note that, still formally, d>H(rt)* = 'l'Jj(rt), and 

5-tj(rt)| =e‘‘"'^+i«'t(r)e-‘‘"'^, «'h(H)L = e‘‘"'^«'(r) (33) 
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3 Operators arising in the GW method for finite systems 

This section aims at providing rigorous mathematical definitions of the operators arising in the 
GW method. For each one of them, we first recall the formal definition given in the physics 
literature, using the second quantization formalism. We then explain how to recast this formal 
definition into a (formally equivalent) satisfactory mathematical definition involving only well- 
defined operators on the fc-particle spaces Hk, with k = 1,N — 1,N,N + 1, the Coulomb space 
C (defined in Section 3.3.1), and its dual C. We finally establish some mathematical properties 
of the operator under consideration, using our definition as a starting point. Unless otherwise 
specified, scalar products and norms are by default considered on Hi = L^(]R^,C). 

We first need to make some assumptions on the physical system under consideration (see 
Section 3.1). We can then define the one-body Green’s functions in Section 3.2. Linear response 
operators are considered in Section 3.3, which culminates with the definition of the dynamically 
screened interaction operator W. We finally introduce the self-energy operator in Section 3.4. 

3.1 Assumptions on the reference A-electron system 

Recall that the reference system with N electrons is described by the Hamiltonian it at on Hn 
defined by (30). Our first assumption concerns the ground state energy of the reference system 
described by it at: 


Hyp. 1: The ground state energy is a simple discrete eigenvalue of itAr. 


In this case, the normalized ground state wave-function of the reference system is unique up 
to a global phase. We also define the energy of the first excited state: 

EIj = min (^a{H]sr)\{E%}y 

Together with we introduce the ground state one-body reduced density-matrix 


Twh,!*') := ^ / 4'^(r,r2,--- ,rAr)«'?,r(r',r2 ,--- ,rAr)dr2---drAr, (34) 

i(R3)N-l 

the ground state density 

Pw(r) = N [ |4'^(r, r2, • • ■ , rAr)P dr2 • • • drAr, 

i(K3)W-l 

and the ground state two-body density 

——[ |4'^(r,r',r3,--- ,rAr)pdr3 ---drAr (35) 

^ J{R3)N-2 

of the reference Welectron system. 

We recall in the following proposition some important properties on 7^, and ^ 
(most of the assertions below are well known; we provide elements of proof in Section 6.9 for the 
less standard statements). Note that both 7^(r, r') and can be seen as the kernels of 

bounded operators on Hi = L^(R^) that we also denote by 7)^, and p% 2 - 

Proposition 26 (Properties of the ground state). Assume that next is of the form 

M 

with Zk S N* and G for all 1 ^ fc ^ M, and that Hyp. 1 is satisfied. Then, 
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(1) the ground state wave-function can he chosen real-valued and S 

(2) the ground state density is in L^(K^,K.) and ^G More¬ 

over, is continuous and everywhere positive on K^; 

(3) the ground state one-body reduced density operator is in 


K-n '■= |7Af G S{T-Li) 0 < 7Ar < 1, Tr?^j(7w) = N, Tr-^i (| V|7Ar|V|) < ooj , 
and satisfies 

V(/,5) G Hi X Hi, (/|7^|ff) = {K\a\g)aim%)n^-, 


(36) 


(4) the kernel 7^(r, r') satisfies the pointwise estimate |7)y(r, r')p < p^{r)p%{r'); 

TV — 1 

(5) the operator p% 2 belongs to S{'Hi), and ||p^_ 2 IIb('Hi) ^ ^—IIPwIIl“- 

Much finer regularity results on are available [13, 14, 48], but are not needed for our 
purpose. Similar results hold true if Vext is replaced by a potential generated by smeared nuclei 
or pseudo-potentials. 

Our second assumption is concerned with the (discrete) convexity of H >->• E^. We assume 
that N ^ 1, and that (with the convention Hq = 0 in the case TV = 1) 


Hyp. 2: E% - E%_, < E%^, - E%. 


In this case, any real number p such that E^ — E^_t^ < p < — E'^ is an admissible chemical 

potential (Fermi level) of the electrons for the ground state of the reference system. The physical 
relevance of this assumption is discussed for instance in [11, Section 4.2]. 

3.2 Green’s functions 

We begin our journey in the GW formalism with Green’s functions. The GW method has been 
designed from the equation of motion for the time-ordered one-body Green’s function G [18], which 
is the concatenation of two meaningful physical objects: the particle Green’s function Gp and the 
hole Green’s function Gh. 

3.2.1 The particle Green’s function Gp 

Rigorous definition of the particle Green’s function. The particle (or forward, or retarded) 
Green’s function is formally defined by (see for instance [12, Section 7]) 

Gpirt, r't') := -i0(t - t') (4'^l4'H(rt)4'][j(rY)l4'^), (37) 

where 0 is the Heaviside function (2), and 4'H(rt) and 'l'Jj(rt) are the Heisenberg representations 
of the annihilation and creation field operators introduced in Section 2.5. As G Hat, we can 
replace 'l'(rt) and '^^rt) by their expressions (33): 

Gpirt, r't') = -i0(t - t') (4'^le“"'^4'(r)e-’(*-*')"'^+i4't(r')e-'‘''f^"l4'^) 

= -i0(t - t') (4'5/l4'(r)e“'*^*“‘')(-^^+^“-®^)4'^(r')l4'^). 

As Gp only depends on the time difference t — t', it is sufficient to study the function Gp(r, r', r) := 
Gpirr, r'O). We then notice that, for all f G Hi, 

[ vI;t(r')|4/^)/(r')dr' = at(/)lvI/^). 
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Introducing 


: 'Hi —?> Hn+ 1 
f ^ aHm%) 

and A+ = (A^)*, we observe that G'p(r, r',T) is formally the kernel of the following one-body 
operator. 

Definition 27 (Particle Green’s function). The particle Green’s function is defined as 

Gp(t) := -i0(r) A;. (38) 

First properties of the particle Green’s function. The study of Gp can be decomposed 
into the study of the operators A+ and The latter is clearly bounded on Hn+i- 

As for the operator A^, we deduce from (32) and (36) that 

{a^f)'S%\a^g)'S%) = (^'^|a(/)a^( 5 )|^'^) = {f\g) - (^'^|a1'(g)a(/)|«'^) = (/|1 - 

or equivalently, 

A+A; = In, - 1%. (39) 

Hence, A^ is a bounded operator from Hi to Hn+i, and A+ is a bounded operator from Hn+i 
to 'Hi. In fact, since 

= (/I(1hi - 7^)1/) = - l%)f\\n, , 

it holds II= 1. The following properties are obtained as a direct corollary of 
Proposition 13. 

Proposition 28 (Properties of the particle Green’s function). The family {Gp{t))t+r defines a 
bounded causal operator on Hi . The real and imaginary parts of its time-Fourier transform are 
in H“®(IR.^,H(Hi)) for all s > 1/2, and are given by 

Re Gp = A+p.v. Im^ =-ttA+H^^+^-^^A;. (40) 

The analytic operator-valued function Gp defined in the upper half-plane by 

y.ev. (41) 

is the Laplace transform of Gp and satisfies 

Gp = lim Gp(.-firy) m H“®(]R;^, H(Hi)) for all s > 1/2. 

77 —>- 0 + 

The imaginary part of Gp is related to the so-called spectral function Ap (see Section 3.2.4). 

Analytic continuation to the complex plane. Let us introduce the particle optical excitation 
set 

Sp := aiHN+i - E%). (42) 

We recall that the operator Hn+i — with domain HAr+inH^(K^(^+^)) is self-adjoint on Hn+i- 
Its essential spectrum is of the form aess{HN+i—E^) = [S 7 V+ 1 , 00 ), and there are possibly infinitely 
many eigenvalues below Sat+i that can only accumulate at S^v+i. According to the HVZ theorem 
[21, 45, 49], T,n+i = E^ — E^ = 0. In particular, Sp is the union of a discrete negative part, and 
the half-line [ 0 ,-boo). 

We next infer from (41) that Gp{z) can be extended to an analytic function from C \ S'p to 
B{Hi). This is of particular interest for the following reason. The operator-valued distribution 
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Gp(w) is highly peaked and irregular (for instance, its imaginary part is a sum of Dirac measures 
on the discrete part of Sp). Instead of studying Gp(w) on the real axis, we will study its analytic 
continuation Gp{z) (defined a priori only in the upper-half plane, but actually on C \ S'p) on the 
imaginary axis p, + iK, where p, < ^ 0 is an admissible chemical potential (see Hyp. 

2). The set Sp can be recovered from uj >->• Gp(/r-|-iw) by locating the singularities of Gp, obtained 
from Gp either by analytic continuation, or by fitting some parameters [36]. We do not address 
this interesting numerical reconstruction problem in the present article. 



Figure 1: Illustration of the analytic continuation: from w i—>■ Gp(a;) to w i—>■ Gp{p -\- iw). 


The following lemma makes precise the behavior of the Green’s function on the vertical axis 
/i -I- iR. It is a direct consequence of the representation (41). 

Lemma 29. Consider p < — E'^. Then the function uj Gp{p + iw) is real analytic from 

to and is in LP(R;^, .S('Hi)) for all p> Moreover, for all w G R, 


Re Gp(/i ict:) — —4l-|_ 


H 


Af-l-l 


-E%- 


uP' + {Hn+1 - E% - pY 


A 


* 

+ 


is a negative, bounded, self-adjoint operator on Hi which enjoys the following symmetry property: 

Vw G R(^, ReGp(/r-I-iw) = ReGp(p — iw). 


For any f G Hi, the function w i—>■ (/|ReGp(^ -I- iw)|/) is non-positive, in L^(R,j), and 



ReGp(^ -I- iw) 


dw 


-7r(/|(lwi -7^)1/)- 


(43) 


The last assertion comes from the spectral theorem, (39), and the equality 


vf; > 0, 



E 

uj^ + E^ 


dw = TT. 


Remark 30. Unfortunately, although ReGp(/r-|-i-) has a sign and (jS) is satisfied for all f G Hi, 


the function uj i—>■ 
fact that 


Re Gp{p -f i’) 


B(Wi) 


does not belong to L^(Rtj). This is essentially due to the 


sup 

E^O 


E 


f;2 


1 

2uj 
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Note that the imaginary part of Gp(/r + iw), 


Im Gp (^ + iw) = —A^ 


UJ 


+ {Hn+1 - E% 


/i)2 


A 


* 

+ ’ 


has no definite sign on and that, for a generic / G 'Hi, the function lo ^ ImGp(/r + iw) / 
does not belong to L^(R^). It will therefore be more convenient in general to work with the real 
part of Gp(ia;) only, especially since the imaginary part can be recovered from the real part (see 
Lemma 31 below). Indeed, the operator-valued functions : uj i—>■ Gp{g, — ij + iuj) are in 
L^(R, S(Hi)) for any ry > 0, and converge to ^ : w !->■ Gp(yi-l-iw) in L2(R,H(Hi)) as 77 —>■ 0+. We 
can therefore apply Titchmarsh’s theorem (see Theorem 9), which gives the following result. 

Lemma 31. Let g, < — E^. The function ^(w) := Gp(/r -|- iw) is the Fourier transform of 

the causal function 

5p(r) = -0(T)A+e-"(^"+i-^"-'^)A;, (44) 

which belongs to L^(Rt-, 5(Hi)). In particular, the Plemelj formulae hold true: 


Re^ = -.^(Im^) and Im^=i 5 (Re^) in L^(R,^,H(Hi)). 


Moreover, the function t i— >■ \\gp{T)\\js{-Hi) exponentially decreasing as |t| — >■ -l-oo. 

Remark 32. The exponential decay of gp is consistent with the analyticity of its Fourier transform. 
This property is of interest when calculating numerically convolutions on the imaginary axis /i-|-iR, 
since convolutions can be replaced, up to a Fourier transform, with point-wise multiplications of 
causal functions which are exponentially decreasing. This approach was advocated in [34], and is 
now routinely used in GW computations. 


3.2.2 The hole (backward) Green’s function Gh 

Definition and first properties of the hole Green’s function. Together with the particle 
Green’s function, we introduce the hole (or backward, or advanced) Green’s function, formally 
defined within the second quantization formalism by 


C/h(H,r't') iQ{t' -t) (4'^|4'j^(r't')4'H(rt)|4'^). 


Observing that 


we introduce 


[ ^{r)\^%)fir)dr = a{m%). 

Jr^ 

A— : Hi Hat—1 

/ ^ a(J)\^%). 


Similarly as before, we note that t/h(rt, r't') only depends on the time difference t — t'. Introducing 
Gh(r, r', r) := GhiiCT, r'O), we see that Gh(r, r', r) is formally the kernel of the following one-body 
operator. 


Definition 33. The hole Green’s function is defined as 

Gh(T) :=i0(-r)Ale‘^(^^-’-'®")A_. 


(45) 


Similarly as in (39), it holds that 

= 7 ^. 

Hence, is a bounded operator from Hi to Hw-i, Af is a bounded operator from Hw-i to Hi, 
and it holds = ||^LHbchw-i.-Hi) ^ 1- The properties of the hole Green’s function 

are quite similar to the properties of the particle Green’s function (compare with Proposition 28). 
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Proposition 34 (Properties of the hole Green’s function). The family (Gh(T))r6R defines a 
bounded anti-causal operator on TLi. The real and imaginary parts of its time-Fourier transform 
are in for all s > 1/2, and are given by 

ReGi; = >llp.v. ( -- -]a_ and Im^ A_. (46) 

V’ - i^N - Hn-i)J 

The analytic operator-valued function Gh defined in the lower half-plane by 

Vz e L, G;(z) := A* \ (47) 

is the Laplace transform of Gh and satisfies 

Gh = lim Gh(-— i?7) in for all s > 1/2. 

77->- 0 + 

Analytic continnation into the complex plane. The hole optical excitation set is defined 
as 

Si,-.= aiE%-HN-i). (48) 

It is clear from (47) that the operator-valued function Gh can be analytically continued to C/S'h- 
Instead of studying the highly irregular distribution ui 1 —>■ Gh(a;), it is more convenient to study 
its analytical continuation Gh on the imaginary axis p. -\- iR, with p > — E^_^. 



Figure 2: Illustration of the analytic continuation: from uj i—>■ Gh(a;) to w 1 —>■ Gh(/r -f iw). 


We can state a result similar to Lemma 29. 

Lemma 35. Consider p > Epf — E'fi_^. Then the function uj i—>■ G\i{p -I- iw) is real analytic from 
Rtj to B{'Hi) and is in LP(R(^, P('Hi)), for all p > 1. Moreover, for all w G R, 


Re Gh (a^ + iw) = Af - 


Hn-1 + P — E^ 






uj^ + {E% - Hn- 1 - pf 

is a positive, bounded, self-adjoint operator, which enjoys the following symmetry property: 

Vw G Rij, ReGh(/.t + iw) = ReGh(Ai — iw). 


For any f G TLi, the function uj ^ if 


|ReGh(At + iw) /) is non-negative, in and 



ReGh(/r -t- icu) 


duj = TT{f\j%\f). 
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The Galitskii-Migdal formula. The hole Green’s function is of particular interest, as it con¬ 
tains useful information on the iV-body ground state. For instance, from the identity = 7 ^, 

we directly obtain Gh(0“) = iy^, so that the expectation value in the ground state of any one-body 
operator (for C S B{'Hi)) can be evaluated via 


K 


N 


2=1 


) = Trn, {Cj%) 


-iTrH,(GGh(0-)). 


This calculation is valid only for one-body operators. It is not possible to obtain the expectation 
value in the ground state of a generic two-body operator from the one-body Green’s function. 
This is however the case for the ground state energy (the expectation value of the two-body 
Hamiltonian Hn in the ground state), as was first shown by Galiskii and Migdal [15]. Alternative 
formulae for the ground state energy are provided by the Luttinger-Ward formula [27] and the 
Klein’s formula [24]. 

Theorem 36 (Galitskii-Migdal formula). For all N ^ 2, the ground state energy can be recovered 
as 


E 


0 _ 

N — 






(49) 

(50) 


The proof of this theorem can be read in Section 6.10. Formula (50) is one way to obtain 
the right-hand side of (49), and is the one found in the original article [15]. There are however 
other ways to obtain (49) from the hole Green’s function, without the use of derivative (which are 
cumbersome to evaluate numerically). One can for instance use the following equality, that we do 
not prove for the sake of brevity, 

Tyhi (^A*_ {Hn-1 + M = Ji^ w^Tr-Hi ^ReGh(At + iw)) . 


3.2.3 The time-ordered Green’s function G 

It is often claimed in the physics literature that the main object of interest is neither the particle 
nor the hole Green’s function, but the function 

g{rt, r't') = t/p(rt, v't') F t/h(ri, r'i'): 

called the time-ordered Green’s function, which can be seen as a convenient way to concatenate 
the information contained in the particle and hole Green’s functions. Obviously, the time-ordered 
Green’s function only depends on the time difference t = t — t' and g{rt, r't') = Qp{r{t — t'), r'O) -I- 
l/h(r(t — t'),r'0). In view of (38) and (45), our definition of the time-ordered Green’s function 
therefore is the following. 

Definition 37 (Green’s function). The (time-ordered) Green’s function is the family of bounded 
operators {G{T))r^s. defined as 

G{t) = Gp(t) -f Gh(r) = -i0(T) A\ + i0(-T) 

The following results straightforwardly follow from Propositions 28 and 34, as well as Lem¬ 
mas 29 and 35. We recall that p, is a chemical potential of the electrons for the ground state of 
the reference system, and that E^ — < fi < E^_p^ — E^. In the following, we introduce some 

G°“(M,j) cut-off functions (j)± satisfying 0 ^ (j)± ^ 1, (j)+ + (p- = 1, Supp((^+) C (E)^ — + 00 ) 

and Supp(0_) C (— 00 , ~ ^n) (®®® Figure 3). These cut-off functions allow us to write prop¬ 

erties of the Green’s function in the time representation without specifying whether t is positive 
or negative. 
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E% - E%_, 


Figure 3: The cut-off functions 4>±. 


Proposition 38 (Properties of the Green’s function). The Fourier transform G = Gp -I- Gh is in 
S('Hi)) for any s > 1/2. The operator-valued analytic function G defined on the physical 
Riemann sheet C \ [Sp U ^h) by 


Vz G C \ (^p U ^h), G{z) ;= A+ 


z - {Hm+i - E%) 




z-{E%- Hm-i) 


A_ 


(51) 


is such that 

lim (/±G(-± iry) = (/±G in i?“®(IR^,S('Hi)) forall s > 1/2. 

The function uj i—>■ G(/i -I- iw) is real analytic from R.;^ to B{TLi), and is in LP(R(^, S(Hi)) for all 
p > 1. Moreover, it satisfies the symmetry property 


Vw G Ri.j, ReG(/r -I- ioj) = ReG(/x — iw). 


For any f G FLi, the function uj i— >■ (/|ReG(^ -I- iw)|/) is in L^(R^), and 



Re G (^ “t" icu) 


dw 


-7r(/|(lHi - 27^)1/). 


3.2.4 The spectral functions ^p, and A 

Spectral functions are essential tools to study many-body effects since they are concentrated on 
(subsets of) the particle and hole excitation sets. 

Definition 39 (Spectral functions). The particle spectral function is the operator-valued Borel 
measure on R^.^ defined by 

V6 G Ap{b) = -ilm^(6) = (52) 

The hole spectral function is similarly defined: 

V6 G Ahib) = hmGi,{b) = A.. 

TT 

The time-ordered spectral function is then obtained as A = Ap -I- Ah- 

With those definitions, the following lemma is straightforward, and is usually referred to as 
the sum-rule for spectral functions (see for instance [11, Section 4.5]). 
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Proposition 40. The spectral functions Ap, Ah and A are S(fHi)-valued Borel measures on 
with supports contained in Sp, Sh and Sp U Sh respectively. For all b S Ap{b), Ah(b) and 

A{b) are bounded positive self-adjoint operators on TLi with norms lower or equal to 1. Moreover, 
0 ^ Ap(bi) ^ Ap{b 2 ) as self-adjoint operators when bi C &2 (and similar inequalities for Ah and 
A), and it holds 

= l-Hi ~ "4h(Rtj) = -4(R(j) = l-Hi • 

Finally, the Plemelj formulae (14) allow us to recover the real part of the Green’s functions from 
the spectral functions: ReGp = 'K^{Ap) and ReGh = •KS){Ah). It therefore holds ReG = tt^A. 


3.3 Linear response operators 

We study in this section the reducible polarizability operator y, which can be defined from the 
so-called charge-fluctuation operator introduced in Section 3.3.1. We give a precise mathematical 
meaning to x Section 3.3.2, and prove Johnson’s sum-rule [23] for x Section 3.3.3. We finally 
define the dynamically screened Coulomb interaction operator (see Section 3.3.4). 


3.3.1 The charge-fluctuation operator pn 

The charge-fluctuation operator is defined, within the second quantization formalism, by (see [11, 
Equation (97)]) 

PB.{vt) := 4'tj(rt)4'H(rt) - p^(r), 
so that the action of this operator on the Wbody ground state is 


PH(rt)lO = (e‘‘(^--0) 

= (4'I(r)4-(r) -p^(r)) 


(53) 


In order to define more rigorously pn, we need to introduce functional spaces of charge densities 
(the Coulomb space) and electrostatic potentials. The complex-valued Coulomb space 


C := {/ G J^'(R3,C) 


/glL(R 3,C), l•l-l/(•)eT2(R^C)}, 


is endowed with the inner product 


(54) 


(/ilDc = 47r / 

Jr: 


/i(k)/ 2 (k) 

IkP 


dk, 


where the normalization condition for the space-Fourier transform is chosen such that its restriction 
to L^(R^,C) is a unitary operator. The space C is a Hilbert space, and L®/^(R^,C) ^ C thanks 
to the Hardy-Littlewood-Sobolev inequality (upon rewriting the products in Fourier space as 
convolutions). The dual of C (taking L^(R^,C) as a pivoting space) is 


C' 


{u G L®(R3,C) 


Vu G (l2(R3,C))^}, 


(55) 


endowed with the inner product 

{Vi\V2)c' ^ [ VTfi-VF 2 = ^/ lkl2^^F2(k)dk. 

47^ Jr3 47r 


We also introduce the Coulomb operator Vc, defined as the multiplication operator by 47rjkj ^ 

1 /2 

in the Fourier representation, and its square root Vc , defined as the multiplication operator by 
(4^)V2|k|-i in the Fourier representation. The following result, whose proof is a straightforward 
consequence of the above definitions, will be repeatedly used throughout this article. 
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1 /2 

Lemma 41. The operator Vc defines a unitary operator from C to C. The operator Vc defines 
a unitary operator from C to Hi, as well as a unitary operator from Hi to C. 

It follows that the adjoint of the unitary operator Vc '■ C ^ C is the unitary operator v* = 
u-i :C' ^C. 

We are now able to reformulate the charge-fluctuation operator in the ground state as a well 
defined bounded operator. For v G it formally holds 


Jr3 / LVi^ / 


vp% 


4'^(ri,...rAr). 


In order to rewrite more rigorously this equality, we introduce the operator 


B: C Hn 

V I—?> 

which is well defined since R) by Proposition 26. In fact, as made clear in Lemma 42 

below, B is bounded. In view of (53), we can finally define the application to 4'^ of the charge- 
fluctuation operator pH(t) as follows: 

PH(t)|^^) (57) 


N 


^u(ri) - {v,p%)c',c 


\K), 


(56) 


Let us conclude this section by giving some properties of the operators introduced above (see 
Section 6.11 for the proof). 

Lemma 42. The operator B defined by (56) is a bounded operator from C to Hn ■ Its adjoint B* 
is a bounded operator from Hn to C which satisfies iI*|'L^) =0. As a consequence, pnl'I'^) G 
L°°(Rt,B{C',HN)), and {ph\K))* e L°°(Rt,B{HN,C')). 


3.3.2 The (symmetrized) reducible polarizability operator y 

Definition of the reducible polarizability operator. The reducible polarizability operator 
x{t, t') is the operator giving the response of the density of the system to perturbations of the 
external potential. It is formally defined by its kernel (see [11, Equation (96)]) 

X{rt,r't') := -i{'i>%\T{pii{rt)pii{r't')}\'i>%)^^ . (58) 

In the above equation, pn is the charge-fluctuation operator whose action on 4'^ is defined by (57), 
and T stands for the bosonic time-ordering operator: 


r{Ai{t)A,{t')} 


Ai{t)A2{fi) ift'<t, 
A2it')Ai{t) if t'> t. 


In view of (57), the definition (58) of the kernel is formally equivalent to the following identity, 
stated for t' < t (a similar equality being true for t' > t): 


R3 




-i / f /(i’)(^w|PH(rt)pH(rY) p(r')drdr' 

-i/ /" /(i')PH(rf) 4 '^dr f g{r')pH{r't')'f>% dr' 

\ JR3 JR3 




Bf 




'Hn 
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In order to interpret x as giving the variation of the ground state density (an element of C) 
generated by a variation of the external potential (an element of C'), we rewrite the scalar product 
in C as a duality braket between C and C: 

{h\f2)c> = {Ti,v-^f2)c,,c- (59) 

This motivates defining x(t, t') as the bounded operator from C to C given by 

In particular, x(t, t') only depends on the time difference t — t', and we write in the sequel x{'^) ■= 
x(t,0): 

x(t) (60) 

It turns out to be useful to symmetrize the action of the polarizability operator using appropri¬ 
ate Coulomb operators. Recall that BvV^ S while G 

Definition 43. The symmetrized reducible polarizability operator Xsym G L°°(Rr,B{'Hi)) is de¬ 
fined by 

Vr G R., Xsym(T) = B*^Bvl/^. 

It is convenient to decompose the symmetrized reducible polarizability operator into two parts, 
namely its causal part and its anti-causal part: 

Xsymir) = xtymM + Xs”ym('r) with x|m('r) = . (61) 

In the above expression, the Hamiltonian can be replaced by 

■“ I {V’ 

This is a consequence of Lemma 42 which shows that Ran (B) C Note that ^ 

~ ^N- 

Properties of the symmetrized reducible polarizability operator. As rigorously stated 
below, the symmetrized polarizability operator has singularities at the energy differences corre¬ 
sponding to excitation energies for a system with a fixed number N of electrons, called neutral 
excitations in [11, Section 8]. We therefore introduce the neutral excitation set 

^0+ := a(i^^^ - E%) \ {0} = a (tJ#, - E%) , 

its reflection Sq := —Sq and Sq ■= SqUSq. Note that C [Ej^ — E^, -boo) so that =0. 

As for Proposition 38, it turns out to be convenient to introduce appopriate cut-off functions. 
Consider (j)\. such that </>)_ and (j)^ are in (^“(Rtj) and satisfy 0 < < 1, -b </>)_ = 1, 

Supp((/)),_) C {-{Elf - A^), -boo) and Supp((/))_) C (-oo, EIj - E%) (see Figure 4). 

Proposition 44. The symmetrized reducible polarizability operator Xsym satisfies the following 
properties: 

(1) (X^ym('r))reR is a bounded eausal operator on Hi while (Xsym(''’))TeK is a bounded anti-causal 
operator on Hi. They satisfy the following symmetry properties: 

Vr G K, Xsym(-r) = Xsym)!-) and xtymM = Xs"ym(-!-); (62) 
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(2) the real and imaginary parts of the time-Fourier transforms of xtyrmXsy-m respectively 
given by 

Reg; . BCA 

and 

In particular, Supp C and Supp C Sq; 

(3) consider the B{'Hi)-valued analytic functions Xsym, X^m and Xs^ respectively defined by 

Vz €C\St, S^{z) := - J . Bvl/^, 

and 

Vz ec\ So, x^{z) := xtymiz) + Xsym{z) = -vf^l'^B* (63) 

It holds 

Vz G C \ Sq , X^in(z) = Xsymi — z) = (^X^in(z)^ 

and 

Vz S C \ Sq, Xsym{z) = Xsym( •?) = (Xsym(-Z)) • 

The functions x^m|u and x^m|L are respectively the Laplace transforms of xtym and x^m; 
and the following convergences hold in S('Hi)) for all s > 1/2; 

liin x4ni(- ± i’l) = X^m, lim (l)±X^{- ± ip) = <('±X^; 

ry-)-0+ 7 )^ 0 + 

(4) for all u! G — E^),Ej^ — E^^, xfff4{aj) = x/^(a;) is a negative bounded self-adjoint 

operator on Hi; 

(5) for all w G K, Xsym(iw) is a negative bounded self-adjoint operator on Hi- 

We omit the proof of Proposition 44 since the first three assertions are similar to those of 
Lemma 38, while the last two ones are direct consequences of (63). 
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On the integrability of As for the Green’s function, uj !->■ is difficult to 

study on the real-axis, and it is more convenient to study its analytical continuation Xsym on 
the imaginary axis iM. This is possible thanks to the existence of the gap {—{Ejf — — 

E^) around 0. The representation provided in Proposition 44 allows one to directly deduce the 
integrability properties of the functions w i—>■ (as in Lemma 29). 

Corollary 45. The functions w i— >■ xi^m(iw) are real-analytic from Ri,, to S{TLi), and are in 
for all p > 1. For any f € Hi, the function uj >—>■ {f\xff^{iu!)\f) is non-positive 
and in L^(R.(^), and it holds 


/ + 00 

(/lA^^(iw)|/)dw = - 27 r 

-OO 


2 

T-Cn 


(64) 


3.3.3 The sum-rule for the reducible polarizability operator x 

The behavior of the reducible polarizability operator in the high imaginary-frequency regime is well 
understood. This asymptotic behavior is given by the so-called Johnson’s sum-rule [23] or /-sum 
rule, the latter terminology being motivated in [11, Section 8.8] by the fact that it can formally 
be seen as some equality involving the first moment of Knowing the large-w behavior of 

is important to design appropriate approximate operators, used in plasmon-pole models to 
avoid the numerical inversion of the dielectric operator (which is computationally expensive). 

The fifth point of Proposition 44 implies that for all uj G the operator —x(ia;) := 

—Vc ^^'^X^^{iuj)vc defines a symmetric, continuous, non-negative sesquilinear form on C': 


V(/,g)GC'xC', 


(/,-x(iw)5)c'.c = 



2{H% - E%) 

iHl,-E%)^+uj^ 



so that, formally, 

hrn (f,-w^xii^)g)c',c = 2{Bf\Hl - El\Bg)n, = 2 (7,vf^B* - E^) Bg) . 

\ \ / /C',C 

The following theorem, whose proof is postponed until Section 6.12, confirms that this limit exists 
and allows one to identify it. 

Theorem 46 (Johnson’s sum rule). The operator 2vf^B*{Hj^ — E^)B is bounded from C to C, 
and 2vf^B*{E[j^ — E^)B = —div (p^V-). Moreover, the following weak convergence holds: 

V(f,g) gC' xC', lim (f,-w^xiiuj)g)^,^ = {J,-div{p%Wg))^,^=f p%Wf-Wg. 

Ul^±CO > > J^3 

For all g G C such that Ag G L^(R^), the following strong convergence holds: 

lim w^x(iw )5 = div (p^Vp) in C. 

0J—¥±OO 


3.3.4 The dynamically screened interaction operator W 

As the name indicates, the two key operators in the GW method are on the one hand, the time- 
ordered Green’s function G, and on the other hand, the so-called dynamically screened interaction 
operator W. The latter operator is defined as 

W{t) = VcSo{t) + vy^Xsym{T)vy^ , (65) 

where Vc is the Coulomb operator introduced in Lemma 41. It is convenient to split W into a 
local-in-time exchange contribution VcSo(t) (although this is not obvious at this stage, (98) below 
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shows that VcSoiT) can be interpreted as an exchange term), and a nonlocal-in-time correlation 
contribution: 


W(t) = VcSo{t) + Weir) with Wc(r) :='UcX(r)wc = (66) 

The properties of the operator Wc(t) S B{C,C') therefore readily follow from the properties of the 
1 /2 

operators Vc and Xsym(T) established in Lemma 41 and Proposition 44. 

3.4 The self-energy operator E 

We give in this section the definition of the self-energy operator S using the Dyson equation 
(see (72) below). Let us emphasize that, while the Dyson equation provides a definition of S in 
terms of Green’s functions, numerical methods work the other way round: an approximation of 
the Green’s function G is obtained from the Dyson equation (72), using an approximation of the 
self-energy operator S. This approach is made precise in Section 4. 

3.4.1 The non-interacting Hamiltonian Hq and associated Green’s function Go 

The self-energy operator is defined as the difference between the inverse of the exact Green’s 
function G and the inverse of some reference Green’s function Go. The reference Green’s function 
is the resolvent of a mean-field non-interacting Hamiltonian. There are several possible choices 
for this operator, discussed in Remark 49 below. In order to remain as general as possible, we 
introduce a one-body operator hi acting on Tti, with domain i7^(]R^), real-valued (in the sense that 
hip is real-valued whenever tp is real-valued), and such that (Tess(hi) = [0,oo). The corresponding 
effective non-interacting Wbody Hamiltonian is dehned on 'Hat by 

N 

Ho,n = y^hi(G)- 

i=l 


We define 


£k ■■ = 


inf sup 

VfcCVfc 


{y\hi\v) 

(u|t;) 


where Vfc is the set of the subspaces of of dimension k. Recall that Sk ^ 0 and that if 

Ek < 0, then hi has at least k negative eigenvalues (counting multiplicities) and Ek is the k^^ 
smallest eigenvalue of hi (still counting multiplicities). We make the following assumption in the 
sequel. 


Hyp. 3: The one-body Hamiltonian hi has at least N negative eigenvalues, and En < En+i- 


This assumption implies that there is a gap between the eigenvalue and the {N+iy*' eigenvalue 
(or the bottom of the essential spectrum if hi has only N non-positive eigenvalues). 

Let us denote by (^i,-- - ,(Pn) an orthonormal family of eigenvectors of hi associated with 
the eigenvalues ei,-- - ,En- Without loss of generality, we can assume that the (pkS are real¬ 
valued. The ground state energy of TLq.at is TIq ^ = ei -I- ... -I- En- The condition en < Eat+i 
ensures that Bq ^ is a non-degenerate eigenvalue of -Bq.n and that the normalized ground state 
= cpi A ■ ■ ■ Arpjv of ffo,N is unique up to a global phase. We introduce the one-body mean-field 
density matrix 

N 

r') := (67) 

This function can be seen as the kernel of the spectral projector l(_oo,/jo)(^i)> where /io is any 
real number in the range {en,£n+i) (it is an admissible Fermi level for the ground state of the 
non-interacting effective Hamiltonian i7o,Af)- The density of the non-interacting system is denoted 
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by pg jy. Results similar to the ones stated in Proposition 26 for pg Jq hold true. Finally, 
similarly as in Section 3.2, we introduce 

and Ao._(/) = a(7)|d>°^). 

Definition 47 (Reference non-interacting Green’s functions). The reference particle, hole and 
time-ordered non-interacting Green’s functions are respectively defined as 

Go.p(t) = -i0(T)Ao,+e-‘"(^«'"+i-<")AS,+ , Go.h(T) = i0(-T)AS,_e‘"(^«'«-i-<«Uo,_, 

and Go(t) = Go,p(t) -h Go,h('r). 

Results similar to Propositions 28, 34 and 38 hold for these operators, but we do not write 
them explicitly for the sake of brevity. However, it should be noted that, in the non-interacting 
case, the Green’s functions have simple explicit expressions in terms of hi (see Section 6.13 for 
the proof). 

Proposition 48. It holds 

Go,p(r) = -i0(r) (l-^, - and Go,h(T) = . 

In particular, for any z € <C\a{hi), 

Go,p( 2;) = (l«i - 7o.w) (2^ - ^i)“^ and Go,hiz) = - hi)~^. (68) 

Hence, 

Go{z) = (z - hi)-^ (69) 

is the resolvent of the one-body operator hi. 

Remark 49 (On the choice of Gq). There are several possible choices for the one-body operator hi, 
although this choice is not really properly discussed in the literature to our knowledge. The first 
option, which is used in the original derivation of the GW method [18], consists in choosing 

hi = ~ 2^ * I ■ I C^O) 

where p% is the exact ground state density. Another option (see for instance [11, page 112[) is 
to consider a one-body operator whose associated ground state density is (as close as possible to) 
the exact ground state density p^. The motivation is that, in this case, the self-energy should be 
smaller. The Kohn-Sham [25] model formally satisfies this requirement. The associated one-body 
operator reads 

hi = ~2^ + P% * I ■ I ^ + Vxc [p^] ) (71) 

where v^c is the (exact) exchange-correlation potential. In practice, approximations of p[^ and 
I’xc [p^] o,fs computed by means of a Kohn-Sham LDA or GGA calculation [25, 32]. This is 
believed to provide a sufficiently accurate approximation of the exact ground state density which 
does not spoil the results subsequently obtained by GW calculations. 

3.4.2 The dynamical Hamiltonian H{z) 

In view of (69), it is natural to introduce the inverse of the time-ordered Green’s function, which 
will correspond to some dynamical one-body Hamiltonian. More precisely, we would like to define, 
at least for each z S C \ K, a one-body operator H{z) such that 

G(z) := — H{z)^ , or equivalently, H{z) = z — ^G(z)^ 

The following proposition, proved in Section 6.14, shows that such a definition makes sense. 
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Proposition 50. Let z G C \ R. The operator G(z) is an invertible operator from 'Hi onto some 
vector subspace D(z) of Hi- Moreover, D(z) is dense in Hi, D{z) C and H{z) is a 

well-defined closed operator with domain D{z). 

Remark 51. We do not know whether the equality D{z) = is true, nor do we know 

whether D[zi) = D{z 2 ) for zi z^- 


3.4.3 Definition of the self-energy operator E from the Dyson eqnation 

We are now able to define the exact self-energy operator E via the Dyson equation. Note that we 
do not define the self-energy in the time domain, but consider only E(z) (as in [11, Section 5.1]). 

Definition 52 (Self-energy). The self-energy operator is defined as 

VzgC\R, E(z) ;=gJ)(z)-1-G(z)-^ = (z-/ii)-( z-H(z)) =H(z)-/ii, (72) 

where hi is the one-body mean-field Hamiltonian introduced in Section S.f.l. 

The operator E(z) is the difference between the one-body dynamical Hamiltonian and the 
reference one-body mean-field Hamiltonian hi. With this writing, E(z) can be seen as the correc¬ 
tion term to be added to the reference one-body Hamiltonian in order to obtain the dynamical 
mean-field one-body Hamiltonian: 

H(z) = fii + E(z). 

4 The GW approximation for finite systems 

4.1 GoW°, self-consistent GW°, self-consistent GW, and all that 

4.1.1 The GW equations 

We now turn to the GW approximation for finite systems. The purpose of the GW approximation 
is to estimate the time-ordered Green’s function G via the Dyson formula (72). Instead of using 
(72) to define the self-energy E(z), we use this equation with some approximation E^'^(z) ofT,{z) 
to obtain an approximation G^'^(z) of the time-ordered Green’s function via 

(g°'^)~'(z) = z- (/ri + EG'^(z)) . (73) 

Using the Dyson equation to define the time-ordered Green’s function is only possible if an alterna¬ 
tive expression of the self-energy operator is available. Such an expression was formally obtained 
by Hedin in 1965 (see [18]). The GW approximation consists in replacing the so-called vertex 
function in Hedin’s equations by a tensor product of Dirac masses. 

The original GW equations were derived on the time domain and on the frequency domain. 
However, as noticed several times in Section 3, the operators involved in the GW equations are 
not smooth on these axes. It turns out that it is formally possible to recast the equations on some 
imaginary axis using Theorem 24. This approach, first introduced by Rojas, Godby and Needs 
[36] (see also [34]), is now known under the name of the “analytic continuation method”. For 
reasons that we will explain throughout this section, these equations are recast as follows within 
our mathematical framework. 
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Definition 53 (GW equations on the imaginary axis of the frequency domain). 
Find + i-) € S('H)) solution of the system 


PGW(i^) = ^^ 1/2 0 GGW(^ + i^') yl/2^ ( 74 ^) 

= (l-H - ^’sym(i^)) -l-Hi, (74b) 

^{ico) = vl/^^iiu;)vl/^ (74c) 

+ = [ G^{^i + iiuJ-uJ')) (D]^^{iuj')duj', (74d) 

J — (30 

G'^'^(/r + iw) = /r + iw — + SGW^^ 0 ^ (74e) 


where hi is the one-body operator defined in (70) and where is the integral operator on Jii 
with kernel 


where 7 q ^ was defined in (67). 

Remark 54. In the GW equations (7)), the chemical potential /i is supposed to be known a priori. 


The GW equations (74) would be the natural equations to work with from a mathematical 
viewpoint (they are formally equivalent to the original Hedin’s GW equations). However, we were 
not able to study (74) for reasons detailed in Remark 55 below. 

As one can directly see, the equations involve quite a large number of operators, which all 
have a physical significance. The operator is the GW approximation of the symmetric 

irreducible polarizability operator, the operator is the GW approximation of the symmetric 
reducible polarizability operator, the operator is the GW approximation of the dynamically 

screened Coulomb interaction operator, and finally jg ^^3 qW approximation of the self¬ 

energy operator. We recognize in Equation (74e) the Dyson equation. The name “GW” comes 
from Equation (74d). 


4.1.2 Different levels of GW approximation 

As mentioned below (see Remark 55), we were not able to study the full self-consistent prob¬ 
lem (74). We will therefore restrict ourselves to the so-called GoW° and GW° approximations. 
We explain in this section how these different models are obtained. 

(i) In the fully self-consistent GW (sc-GW) approximation, we assume that the full problem (74) 

is well-posed, so that there exists a (unique) solution G®^. It is then solved self-consistently: the 
idea is to start from some trial Green’s function, and keep updating it with (74) until convergence. 
This method is for instance used in [7, 8 , 26, 37, 41]. It was implemented only quite recently due 
to its high numerical cost (one needs to perform the inversion in (74b) at each iteration). 

(ii) In the so-called self-consistent GW° approximation, or simply GW° approximation, only 
the Green’s function (and not the screened Coulomb operator) is updated in (74d) (see for in¬ 
stance [42, 46]). This partial update not only speeds up the calculation (the inversion in (74b) is 
only performed once), but is sometimes in better agreement with experimental results than the 
sc-GW approximation. This is the model that we study in Section 4.3. 

(iii) Finally, most works simply consider the GqW^ approximation, where only one iteration 
of the sc-GW (or equivalently one iteration of GW°) is performed. This model is very popular 
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due to its relatively low computational cost, and provides already very satisfactory results (see for 
instance [4]). 

Let us also emphasize that it is unclear that a solution of the fully self-consistent GW model 
is a better approximation in any sense to the exact Green’s function than a non self-consistent 
approximation such as the one obtained by the GqW*^ approximation. This is discussed in [11, 
Section 9.8], where the author also comments on the possibilities to update the effective one-body 
operator hi + along the iterations. 


Remark 55. We do not know how to give a proper mathematical meaning to Equation (Ha). 
More specifically, one would like to define, for a reasonable choice of Green’s function G’^pp, the 
operator 

Vw G ]R„, Psym[G'’'PP](iw) := (^J G^PP{p + i(tLi -|- w')) 0 G’‘pp(^ -I- iw') dw'^ 

and we would like this operator to be a self-adjoint bounded negative operator on T-Li. It is the case 
for instance when G^-pp is the non-interacting Hamiltonian Go defined in (69) (see Proposition 59 
and Remark 64), or when G'^pp is the exact Green’s function defined in (51) (this fact can be proved 
by adapting the arguments given in Section 4-2.2). We were not able to obtain this result for a 
generic class of approximate Green’s functions G^-pp, say G^pp of the form (74e) with 
in a small ball of L°°(M.cj,B{'Hi)). 

For this reason, we will not study the self-consistent GW equation (74). 

4.2 The operator and the random phase approximation 

The remainder of this section is devoted to the study of the GW° approximation (which includes 
the GoW° approximation), which amounts to study the two equations (74d)-(74e) with a specific 
fixed choice of the screening operator . This approximation bypasses the difficulties mentioned 
in Remark 55. In order to present and study the GW° approximation, one must first define the 
operator W^. 

4.2.1 The RPA irreducible polarizability operator 

The GW approximation of the irreducible polarizability operator P is formally defined as 

P®'^(r, r', t) = -iG(r, r', r)G(r', r, -t). (75) 

When the Green’s function G is the non-interacting one Gq defined in (47), this also corresponds 
to the so-called random phase approximation of the reducible polarizability operator (compare for 
instance (83) with the expression in [6]). We therefore define 

P°(r, r', t) := -iGo(r, r', r)Go(r', r, -t). 

This operator is expected to have properties similar to the operator y dehned in Section 3.3.2. In 
particular, P^{t) is expected to be a bounded operator from C' to C. It is therefore more convenient 
to work with its symmetrized counterpart Psyni(T) := which is expected to be a 

bounded operator on Hi. It is possible to decompose Pgym -^sym = -^sym + -^sym where, using 
the kernel-product 0 defined in Section 2.4, and the explicit expressions of Go,p and Go,h given in 
Proposition 48, 

iOiryj^GoAr) © GoA-t^J^ 
i0(T)wy2 ((Iwi - iIn) 0 To.ive'^'*') 
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(76) 

(77) 


and 


^’sym('r) = -i0(-T)vy^Go,h(T) © 

= 0 - 7o%) ■ 

Actually, with this definition, we were not able to give a meaning to Pgy^ (it may not be a bounded 
operator on'Hi). We therefore prefer to use the modified kernel-product 0 defined in Remark 19. 
Our correct mathematical definition for Pgy’” then is 

Psy^ir) = -i0(-r)uy2Go,h(r) 0 Go.p(-r)uy2 (78) 

= -i0(-'r)t'y^ (7aAre“'^^'0 {tui - 7 o,n) (^9) 

As will be shown in Lemma 56, this amounts to defining P°’“(r) = P°’+(—r). We recall that 
7 g ^ is the orthogonal projector on the vector space spanned by the eigenvectors of hi associated 
with the lowest N eigenvalues (see (67)), so that 

N 

7o.iv = XI 

fc=i 

where hi(j)k = £k4‘k, and the eigenfunctions (j)k are real-valued and orthonormal. The following 
result shows that our definitions make sense, and gives explicit formulae for P*^’+ (see Section 6.15 
for the proof). 

Lemma 56. The family (Rsym ('^))reR, de/ined by (76) is a hounded causal operator on 'Hi, 
while (T’sym(^))TGR ('^^) ® bounded anti-causal operator on 'Hi- It holds Pgy’“(T) = 

■P°ym(-'r) and 


N 

Psfmi'r) = -i 0 (T) X (l-Hi - tIn) {in, - 7S.iv) 




(81) 


Remark 57. For 1 ^ fc ^ N, the notation (fk in (81) refers to the multiplication operator by 
the function 4>k- It is a bounded operator from C to 'Hi, and from 'Hi to C (see the proof of 
Lemma 56). The operator is bounded on 'Hi, and one can check that its adjoint on 'Hi 

is {(fkvV^)* := 

The properties of the Laplace and Fourier transforms of Rgym easily deduced from (81) 
using Proposition 13 and Lemma 14. 

Proposition 58. The function z i—> Psy’m)^:) is analytic on the upper half-plane U, and can he 
analytically continued to the lower half-plane L through the semi-real line {—oo,eN+i — £n)- For 
all z € C \ [ffTv+i — £n, oo), 


N 


Ps°yi(^)=X^X<^fe ( ^ 


l-Hi 


7o,n 


k=l 


z — hi -\- Ek 


4>kv, 


1/2 


(82) 


Moreover Psym(- + ii?) converges to Piym 


in H ^(R,^,B{'Hi)) as p ^ 0+, with 


N 


Re Psym = P.V. 




Iwi — 


^0 

10,N 


\k=l 


— hi -\- Ek 


(fkV, 


1/2 


and 


N 


Im R 


00 




\k^\ 
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It also holds 


Vz e C \ [eat+i - Eat, Oo), Psym{z) = Psym{-z), 

SO that, for z S U U L U (—(eat+i — Eat), Eat+i — £n), 

_ __ N / h - \ 

Psymiz) = “ tS.Tv) j “ 7o.iv)</>fe^y^- (83) 

The properties of Psym and of Psym can be directly read off from the previous expressions. For 

instance, we see that ImPsym and ImPsym are negative operator-valued measures, with support in 
(ea/+i —Eat, oo) and (—00, —(ea/+i —Eat)) respectively. For w in the real gap (— (eat+i — Eat), Eat+i — 

Eat), we see that Psym(w) = RePsym(w) is a negative bounded self-adjoint on Hi. 

For our purpose, we only need to know the behavior of Psym on the imaginary axis iK^^. We 
summarize the corresponding most important results in the following proposition (see Section 6.16 
for the proof). 

Proposition 59. It holds 


— ^ / 
Vw G Ps°m(iw) = -2^ (l^i - 7o,Ar) ( 


hi - £k 


oP -b (hi - Ekf 


(l-Hi -7o,Ar) 


_ (84) 

In particular, for all to G the operator Pg‘^m(i‘^) a negative, self-adjoint bounded operator 

on Hi satisfying Pg'^m)”!*^) = -Psym(i'^)- addition, the function uj 1 —Pg‘^m(ia^) analytic 
from R^ to 5 (^ 1 ), and is in , S(Hi)) for all p > 1. For any f G Hi, the function uj 1 —>■ 


■Psvm(iw) 


fj is non-positive, in L 

P+OO 




), and 

/) dw = -27r (/ vl/^ ((1-Hi - tS.at) 0 7o,n) 

N 

= -2n{fJ2 (l«i - 7o%) <l^kvy^ 


k=l 


Finally, there exists a constant C G R'*' such that 


Muj G 


0 < -^s°ym(i‘^) ^ 


c 


(W 2 + 1 ) 


1/2 ( 


,,l/2„0 1/2 

Po,N^c 


). 


(85) 


( 86 ) 


where pg jy is the multiplication operator by the (real-valued) function PQpf, the latter operator 
being bounded from C to C. 

The sum-rule for the operator P°. We end this section with the sum-rule for the operator 
po = which goes from C to C. We postpone the proof until Section 6.17. 

Theorem 60. The operator ~ 7o n)(^i ~ £k)4>k Is bounded from C to C, and it 

holds 

N 

2 ^kit-Hi - 7o.Ar)(^l - £fe)<(>fe = div (po,ArV-)- 

k^l 

Moreover, the following weak-convergence holds: 

'i{f,g)GC'xC, hm (j-ujWoiiuj)g) = (/,-div ^ • Vp. 

cl)->-±oo \ / C',C ' / c 5 

Finally, for all g € C such that Ag G L^(R^), the following strong convergence holds: 

lim w^P°(ia;)(? = div (pg ^Vp) inC. 

CJ—>-±oo ^ ' 
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This sum-rule automatically leads to a sum-rule for the reducible polarizability operator in the 
random phase approximation (see Theorem 67). 


4.2.2 The analytical continuation method 

In this section, we explain why (74a) can be thought of as a natural reformulation of the usual 
physical definition (75), and why problems arise with Definition (74a) (see Problem 55). This 
section also serves as a guideline to understand why (74d) is a natural reformulation of the usual 
physical definition of (see (97) below). In the previous section, we gave the properties of 
using the explicit expression of given in (81). While this approach simplihes the proofs, it 
somehow hides some structural properties that we highlight in this section. 

Recall that = pOy’+ -h with 


and 


where 


-Psym(T) = -i0(T)^^c''^G'o,p(T) 0 Go.h(-T-)uy2 


^sym(T) = -i0(-'r)i^y^Go,h(r) 0 Go,p(-r)uy^ 


Go.p(r) = -i0(T)7lo.+e-‘"(^‘>'"+i-<«)^S.+> Go,h(T) = i0(-T)7l*__e'"(^»'«-i-<«Uo._. 

The idea is to use the results of Theorem 24. We first consider P^^^, and prove that the hypotheses 
of Theorem 24 are satisfied. This is given by the following lemma. 

Lemma 61. There exists a constant C G K’*' such that, for any f G TLi, it holds Aq,- G 

62 ( 7 ^ 1 ) with 


An- 


- {'-'Tf) 


BaCHi) 




Moreover, Hq^n+i - E° j,j > en+i and Hq^n-i - pQ.iv ^ -^n- 

Proof. The first point comes from the fact that Aq _Aq^- — 7o n tha,t f € C’ ^ 
whenever / G TLi, together with Lemma 77. □ 

In particular, the hypotheses of Theorem 24 are satisfied, and we deduce that for any v' > ejv 
and V + v' < £n+i, 


Vw G 


p0,+ 

JA sym 


1 p +00 __ ___ 

(zz -I- iw) = — J ^Go,p(i^ Pv' P i(a; -I- w')) 0 Go,h(i^' + iw')J . 


(87) 


We treat Psym is a similar way, and hnd that for any v' < e^v-i-i and v P v' > £j<s, 


Vw G 


P 


0 ,- 

sym 


1 /■+°° /__ ___ \ 

{v p\uj) = — J v]J'^ (^Go,h(j^ -I- zz' -I- i(a; -I- J f) 0 Go,p(i^' P ia;')J . 


( 88 ) 

Actually, the kernel-product 0 in the latter expression can be transformed into the kernel- 
product 0, thanks to the following lemma, whose proof is given in Section 6.18. 

Lemma 62. For any v' < £n+i, any v P v' > En and any lo,uj' G 

Go,h(^' P E P i(w -|- w^)) 0 Go,p(z^^ P itu'') = Go,h(^^ -I- zz^ -I- i(w -I- lu')') 0 Gojp(z^^ -I- iw^), 

as bounded operators from C to C. 
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We can perform the same type of calculation for Gh 0 Gh- Following the proof of Theorem 24, 
we deduce from G\i{t) 0 Gh(—r) = 0 that, for any v' > en and v + v' > En, 

1 /•+°° /__ __ \ 

Vw e R„, — j \Gq^y,{v + v'+ i(uj + uj')) QGq^\,{v'+ = Q. (89) 

Similarly, from Gp(t) 0 Gp(—r) = 0, we deduce that, at least formally^ for any v' < En+i, and 
any v + v' < En+i, 

1 r-\-oo \ 

Vw e ]R„, — j ^Go,p(z/+ ^'+ i(a; + w')) 0 Go,p(i^'+ ia;')J uy^di:^'= 0. (90) 


Remark 63. The last equality is formal, in the sense that the integrand Go,p 0 Go,p is actually 
not well-defined: it does not define a bounded operator from C to C. However, we can proceed as 

follows. For u) G R^j, let ^^’’^(ia;) be the operator defined on the core TLi nC by 


yf,gGninc, (/ 




iw) 


r*+oo 




Co.pG + r + i(w + w')) (vy^g'j Gop(E + iw') (^y^/) 


dw'. 


Noticing that vl^'^ f and vV^g are in Hi since f,gGC, and reasoning as in the proof of Lemma 62, 
we can prove that the operator in the trace is indeed trace-class, with 


Tr-Hi Go^p{i' + i/' + i(a; + uj')) (vl^'^g'j Go,p(j^' + iw') 


where p^ is an integrable function independent of f and g. Moreover, following the proof of 
Theorem 24, we can prove that, as expected. 


v/,gGHinc, 



= 0 . 


The unique continuation on Hi of therefore is the null operator. It is unclear to us how 

to extend a similar reasoning for a generic class of approximated Green’s function G^-pp. 

By gathering (87), (88), (89) and (90), we find that, for any u' G (£:Ar,£w+i) and n -\- N G 

{en, ew+i), 

Vw G Rt,;, P°y^{iy-\-iuj) = — J ^Go(i^ + + i(cLi'+ w)) 0 Go(i/'+ ia;')j uy^dw'. 

In particular, this equality holds for the particular choice v' = po and ly = 0. 

Remark 64. To summarize the work performed in this section, we transformed the equation 

P^{r, r', t) := -iGo(r, r', r)Go(r', r, -r) (91) 

into: for any v' G [en,£n+i) ^ £ {^N — i^',£n+i ~ ^') 

P'^{v-G\-) = — J ^Go(j^ + + i(a;'+ •)) 0 Go(!^'+ iw')j dw'. (92) 

Note that the manipulations performed in this section to transform (91) into (92) are possible 
since the two operators involved in the kernel-product (here, both are equal to G^{z)) are analytic 
on some common domain U U L U (a, 6 ) with a < b (the presence of a gap is important to deform 
the contour as in Theorem 24). 


38 












4.2.3 The RPA reducible polarizability operator 

In order to calculate the GW approximation of the self-energy, one needs the reducible polariz¬ 
ability operator y, defined in Section 3.3.2. Unfortunatly, the expression of y is not accessible 
in practice. One needs to approximate this operator. The GW approximation, which amounts 
to approximating the so-called vertex function, provides a natural approximation y®'^ of y: in 
Equation (74b), y^'^ is defined from (see also [11, Equation (103)] or [18, Equations (A.20) 
and (A.28)]). However, in view of Remark 55, the definition of y®'^ is not well-understood math¬ 
ematically. In the GW° framework, we use the RPA reducible polarizability operator y°, which is 
itself defined in terms of the RPA irreducible polarizability The GW° approximation of the 
(symmetrized) reducible polarizability operator is usually defined in the frequency domain as 



The formal analytic continuation of the above definitions is (see [11- Equation (139)]) 

^(z):=(l„,-J^(z)) (93) 

Note that we use the “tilde” notation in ygy^n, although it is unclear that this operator-valued 
function is indeed the Laplace transform of some operator-valued function in the time domain. 

Also, it is a priori unclear whether the operators 1 -Hi — fsym {^) oi' ~ T^sym invertible. 

This is however the case for appropriate values of z, as shown by the following lemma. 

Lemma 65. For z G {—{sn+i — eN),£N+i — £n) and z G iK, the operator — Pgy^{z) is 
invertible. 

This result is a direct consequence of the explicit formula (83) for which ensures that 
'fsym('^) i® ^ bounded self-adjoint negative operator for the values of 2 under consideration. Let 
us deduce some extra properties of y^. 

Lemma 66. For any w G M, the operator yOy,jj(ia;) is a bounded, negative, self-adjoint operator 
on Hi, satisfying y0yjj,(—iw) = yOyj^(iLti), and such that 

Psymi^oj) ^ y°ym(iw) < 0. (94) 

The function uj i—>■ y0y,j,(ia;) is analytic from to S{Hi) and is in , S{Hi)) for all p > 1. 

Finally, there exists a constant C G K.’*" such that 

0 ^ -X°ym(iw) ^ (^2 -^1)1/2 Po,Nnl^^) ■ (95) 

This result is deduced from the definition (93), the inequality a: ^ (1 — x)~^ — 1 ^ 0 for a; ^ 0, 
and Proposition 59. 

Sum-rule for y°. From the sum-rule stated in Theorem 60, we readily deduce the sum-rule 
for y*^ := Vc ^^^Xsym^c which is a bounded operator from C to C. Indeed, from the equality 
(1 — x)~^ — 1 = X x^{l — a;)“^, we obtain 

VwGK,^, yOy,„(ia;) = P,%(ia;)-f (P,%(ia;)) - P0y,„(ia;)) . 

In particular, 

Vcc G M., ccV(icc) = cc2p0(ic^) + ^ (c^2p0(icc)) _ i^(icc))”' (co^P^iiw)) . 
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This shows that the asymptotic behavior of x°(iw) is, at dominant order, the same as for P0(ia;). 
Taking the limit uj —>■ ±oo leads to a theorem similar to Theorem 60, whose proof is skipped here 
for the sake of brevity. 


Theorem 67. The following weak-convergence holds: 

'i{f,g) eC xC', lim //,= (7, -div ^Vg)) = 

CJ—\ / C' ,C ^ / L. ,t- 


/ 

Jr3 


Po.ivV/ • Vg. 


For all g € C such that Ag S 


the following strong convergence holds: 


lim a;^x°(iw)5 = div (po in C- 

CJ— >-±oo ’ 


By comparing Theorems 67 and 46, we see why using (71) instead of (70) for the definition of 
hi may lead to better approximations, since Pq]^ = p% in this case, so that the GW approximation 
xGW ^ becomes exact in the high imaginary frequency domain. 


Theorem 67 is useful for the design of the so-called Plasmon-Pole models (PPM) [22, 47, 16, 
10]. Since the definition (93) requires the computation of a resolvent, the calculation of x^i^) is 
numerically very expensive in practice. Some authors suggested to approximate x° by an operator 
^PPM -^hich is computationally less expensive. In practice, x^^'^ bas a prescribed functional form, 
with adjustable parameters. Different approaches are taken in order to tune these parameters, and 
the previous sum-rule provides a standard way to fit some of them. This is done for instance in the 
PPM by Hybersten and Louie [22] and in the PPM by Engel and Farid [10]. In the later article, 
the authors extensively comment on the fact that this sum-rule is an important requirement to be 
satisfied for a PPM. 


4.2.4 The RPA dynamically screened operator 

From the approximation x'^ of x, we directly deduce the approximation of W. Following the 
path taken in Section 3.3.4, we define 

W°{z) ■.= Vc + \^iz) with wj^iz) (96) 

This operator, when well-defined (say on the gap (—(ejv+i — Eat), Eat+i — Eat) or on the imaginary 
axis iK) is a bounded operator from C to C'. The properties of are directly deduced from the 
ones of Xsymi we do not repeat them here for brevity. 

4.3 A mathematical study of the GW° approximation 

4.3.1 The GoW° approximation of the self-energy 

In this section, we study the GoW° approximation as a preliminary step to the study of the self- 
consistent GW° approximation. This will help us understand some technical points to address in 
the analysis of the GW° method. 

The GqW® approximation of the self-energy operator is formally defined as 

E°°(r, r', t) := iGo(r, r', r)IT°(r, r', -r+). (97) 

Here, Gq represents the Green’s function of the non-interacting system introduced in Definition 47, 
and is the random phase approximation of the dynamically screened operator defined in Sec¬ 
tion 4.2.4. Already one difficulty arises: in Section 4.2.4, we only defined the function W^(z) 
on the complex frequency domain, but we did not define some operator-valued function on the 
time-domain. In this section, we assume that the function W^(z) is indeed the Laplace transform 
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of some operator This will allow us to transform (97) into a formally equivalent definition 

that only involves W'^. The resulting definition will be our starting point for the GW° approxi¬ 
mation. 


With the kernel-product defined in Section 2.4.1, the definition (97) can be recast as 

S°°(t) =iGo(r-)©W^°(-T). 

In view of the decomposition provided in (96), it is natural to split into an exchange part E™ 
and a correlation part S™ (the terminology is motivated below): 

= E°° -h E™ with E°°(r) = iGo.h(O-) © Vc5o{t) and E™(t) = iGo(r) © IFc(-r). 

Let us first consider the exchange part. As iGo,h(0“) = —Joiy, we obtain 


E^'^(T)=if,5o(T), 

where is the integral operator on Hi with kernel 

1 7S.iv(r,r') 

K^{r,r) := --—. 


(98) 


(99) 


We recover the usual Fock exchange operator associated with Jq which justifies the terminology 
“exchange part” for E°°. Let us now consider the correlation part. Observing that 

• Go is analytic on U U L U (sjv, Sn+i) (hence has a gap around /j,o) ; 


• is analytic on U U L U (—(eat+i — Eat), Eat-i-i — sn) (hence has a gap around 0), 

we can use the same ideas as in Section 4.2.2. By analogy with Remark 64, we recast the physical 
definition of E*^*^ in (97) in a formally equivalent definition in the complex frequency plane. This 
reformulation was first given by Rojas, Godby and Needs [36] (see also [34]), and is now known as 
the “contour deformation” technique. 

Definition 68 (GqW^ approximation of the self-energy). The exchange part of the self-energy in 
the Go approximation is defined in the complex frequency domain by 

Vz G C, E™(z) = 


while the correlation part is defined, for v' G (—(eat+i — Eat), Eat+i — Eat) and v -\-v' € (eat, Eat-i-i) 
by 

- 1 ~ - 

Vw G Ktj, E[l°(i/-I-iw) = — — / Go {vi{ujuj')) QW^{r'ioj') duj. 

J — OO 

The fact that the above quantity is independent of the choice of v' comes from the analyticity 
of the integrand on the region of interest. In practice, we will focus on the case i/' = 0 and r = po^ 
and therefore consider the function 9 w i—>■ Eg°(/ro + iw) defined by 

Vw G Rt,;, EOO(/xo-I-iw) =-— / Go (/.to + i(w-I-w')) © ^^(ia;') dw. (100) 

The next proposition shows that the above definition makes sense. 

Proposition 69. The operator arising in the exchange part E™ of the self-energy is a negative 
Hilbert-Schmidt operator on Hi- Furthermore, for any w G R^,,, the operator E™(/to -I- iw) is a 
bounded operator onHi, and satisfies E™(/to —iw) = E™(/io+iw)*. The function tv i—>■ E™(/to-|-iw) 
is analytic from R(^ to and is in L^’(R^, S(Hi)) for all p > 1. 

The first statements of Proposition 69 can be seen as a special case of Proposition 73, while 
the symmetry property for the adjoint and the integrability follow from the properties of G*^ 
and W^. 
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4.3.2 Well-posedness of the GW° approximation in the perturbative regime. 

We finally study the GW° approximation. Following our definition (100) of the GoW° approxi¬ 
mation of the self-energy, we recast the GW° equation as follows. 

Definition 70 (The GW° problem on the imaginary axis in the frequency domain). 

Find G^'^° S L°“(R(j, 0 ( 7 ^ 1 )) solution to the system 




+ 


r*+oo 


= K^- 


+ i(a; + w')) 0 WO{ioj') dw', 


GGW“(^o+i^) 




+ iuj — i hi + SGwo 


1 -1 


where /ii is the one-body mean-field Hamiltonian defined in (70) and is the exchanqe operator 
defined by (67)-(99). 

Remark 71. We are looking for a solution in Note that the true Green’s func¬ 
tion G(/i -I- i-) is in BfiHi)) for all p > 1 (in particular for p = oo). We chose to work with 

L°°(]&i^,B{'Hi) for simplicity, but it is possible to work with other spaces LP(Rtj, 5 ( 7 ^ 1 )) withp > 1. 


Since this problem seems quite difficult to study mathematically, we will only study it in a 
perturbative regime. More specifically, seeing as a correction term (see the discussion after 

Definition 52), we propose to study the following problem. 

Definition 72 (The GW° problem on the imaginary axis on the frequency domain). 

Find S B{'Hi)) solution of the system 


(GWl) I 


( -- 1 ^- 1-00 

\ ^^^°ipo+iuj) = K,-— i 


GGW“.(^o + i^) = 


G®^° {po -\- i(w -I- ui')) 0 VF°(iw') dw' 

n-l (101) 

GW« / 


^0 + iw - [hi -\- ^ {po -\- iw) 


According to (69), the unique solution for A = 0 is the Green’s function for the non interacting 

system G^^°=o = Gq. This fact will allow us to treat the equation perturbatively. The exact GW° 
equations correspond to the case A = 1. Of course, several other choices of perturbation can be 
used. For instance, we can put the parameter A in front of the correlation part of the self-energy 
only. This amounts to considering the Hartree-Fock Hamiltonian as the reference Hamiltonian 
(instead of the Hartree Hamiltonian). The theory that we develop here can be straightforwardly 
generalized to such other cases. 


ft is convenient for the mathematical analysis to introduce the functionals s and g respectively 
defined as 


s: L^R^,B{ni)) ^ L°°(R„,H(77i)) 


«+oo 


G^PP{po + i-) hO s G^PP {po + i-) := K,, - — G^PP{po + i{-+ uj')) Q W;){iuj') duj' 
L J ZTT 


and 


Ba: L^(Ru.,BiHi)) ^ L\R^,B{ni)) 

lW^{po+i-) ^ (^o + i.) := J^o + i._(^/ii+AE55(^o + i.)) 


With this notation, G®'^" is a solution of the GW° equations (101) if and only if it is a fixed-point 
of Ba o S- The fact that these maps are indeed well-defined is proved in the following proposition 
(see Section 6.19 for the proof). 
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Proposition 73. The operators is a bounded linear operator from to 

On the other hand, for all M > 0, there exists Am > 0 and Cm G K.'*' such that for all 0 ^ X < Am, 


and all E^-pp such that E^-PPi/io + i-) 

L“(R„.Z 3 (Wi)) 

defined as an element of BifHi)) (1 L°°(Ri^,B{'Hi)), with 


0A 


l]app 


ih-o + i-) 


0A 


< M, the function (//q + i’) is well- 

), with 

E'^PP ifio + i-) _ ^ Cm- 


Moreover, for all G L°°{R^,B{'Hi)) such that + i’) 

1 ^ j ^ 2 , 

0A 


.M-Hi)) 


^app 

-0A 

^app 

II 

^app 

- EfP) gA 

^app 


^ M for 


( 102 ) 


To prove the existence of a fixed-point for Oaos, we rely on Picard’s fixed-point theorem. Since 
the solution of the GW°^q equations ( 101 ) for A = 0 is Cq, we are lead to introduce, for r > 0 , 
the (closed) ball 




^Go, r^ — ■|G’^PP(/ro -|- i-) G L^(R,j, BiTLi)), G®'PP(/io -I- i-) — Go(/io + i 


^ r 


The existence of a fixed-point is given by the following theorem (see Section 6.20 for the proof). 
Theorem 74. There exists A* > 0 and r > 0 such that, for o/Z 0 ^ A ^ A*, there exists a unique 
element G S ^Co,r^ solution to the GW° equations (101), or equivalently to the fixed point 

equation 

^1=qxos(g^I^ . 


In addition, for all uj G -I- iw) is an invertible operator, and 


^GGWx (^g-p _^Go(p(o + i-) 


-1 


< oo. 


(103) 


Finally, the iterative sequence (gA o s)^ 
and C G K'*' such that 


GGWO_(gAos)* 


Gr 


Gn 


Z,”(R,.,.B(Wi)) 

converges to G^'^°>-, and there exists 0 ^ a < 1 


(Mo + i-) 


^ CaC 


..BCHi)) 


Remark 75. It is not difficult to deduce from (103) that the function ui i—>■ GG''’”(^ -p ioj) is 
actually in LP(R,^,B{'Hi)), for all p > 1. 


5 Conclusion 

This article is, to our knowledge, the first attempt to formalize with full mathematical rigor the 
GW theory for finite molecular systems derived by Hedin in his seminal work published in 1965 [18] . 
In Section 3, we have provided a mathematical definition of some one-body operators arriving in 
many-body perturbation theory for electronic systems, namely the one-body Green’s function G, 
the spectral function A, the reducible polarizability operator y, the dynamically screened interac¬ 
tion operator W, and the self-energy operator E. 

In Section 4, we have worked out a mathematically consistent formulation of the GW° ap¬ 
proximation of the GW equations, and we have proved that the GW° model has a solution in a 
perturbation regime. As a by-product, we have also shown that the widely used GoW° approxi¬ 
mation of the self-energy makes perfect mathematical sense. 
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6 


Proofs 


6.1 Proof of Lemma 3 

Let s > 1/2. For / S and ip S 


{TTf,p)^, 


^',5^ 


= 27r 


/(-t) 


= 27r 


f{-T){FT^p){T) dr 


( 1 + t 2)'>/2 


(l + r2)*/2(J--V)(r)dr 




H’ 


where we have used the Cauchy-Schwarz inequality in the last step. By density, J-xf can be 
extended to a linear form on iL®(R). The equality case WJ-tIWh-” = C'sH/IU” is obtained for 
constant functions. 


6.2 Proof of Theorem 10 

Proof of (i). The analyticity oig directly follows from the results of [38, Chapter VIII]. 


Proof of (ii). Let s > 1/2, and consider p G Relying on the fact that g{- + if]) can be 

seen as the Fourier transform of r i—>■ g(T)e~^'^, we obtain 


{gi- + ^V), v)h-pH‘' - (g, <f)H-pH‘ 


{ge-^\p} 


f 


gi^) 


(l + r2)V2 


) (1 + r2)®/2^(r) 



l) dr. 


(104) 


where the integral makes sense since r i—>■ ( 7 (r)(l + r^)“®/^ and r i—>■ {\ + T‘^)''^/'^ip(T) are in L^(]R). It 
is then possible to extend the above formula to any p G Moreover, by the Cauchy-Schwarz 

inequality, 

\{g{- + iv),p)H-pH‘' - {g,<f)H-pH- \ ^ Irj,s\\p\\H‘\\g\\L^, 


where 



( l - e -''^)2 
(1 -I- r2)s 



1/2 


< oo. 


Therefore, || 5 (-+ if?) — ^ ||<?|lL“lr/,s- By dominated convergence, Irj,s —>■ 0 as ?? —>■ O'*', 

which allows us to conclude to the strong convergence ofg{- + if?) to 'g in iL“®(IR.^). 

A similar computation shows that, for 0 < ryi ^ f ?2 and s G K., 


||5(-+if?i)-ff(- + if?2)||ff. 





1/2 


where we crucially use that 7 ?i > 0 to ensure the convergence of the time integral when s > — 1 / 2 . 
The right-hand side goes to 0 as f ?2 goes to f?i by dominated convergence. This allows one to 
conclude to the continuity of i? i—>■ g{- + if?) from (0, -|-oo) to iJ'*)®). When s < —1/2, it is possible 
to pass to the limit f?i —>■ 0 and obtain the uniform continuity from [ 0 , -l-oo) to 7J'’(R). 


Proof of (iii). We follow the approach used in [43] for instance. Fix zq G U, and consider, for 
R > 0 and 0 < a ^ Im(zo)/2, the oriented contour ^ in the complex plane composed of the 
semi-circle ia -I- i?e‘® for 0 ^ 0 ^ tt and the line ia + uj for —R ^ uf ^ R. The value R is taken 
sufficiently large for zq to be inside the domain encircled by the contour (see Figure 5). 

By Cauchy’s residue theorem. 


g{zo) = 


2m 


Z — Zq ZITT 


1 

ifT J-R 


gjuj + ia) 
j^u; + ia — Zq 


dw 


1 

27r 


g{ia + Re'^) 


R&' 


i?e'® -I- ia — zo 


d 6 ». 

(105) 
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Now, for z G U, 


so that 


\9iz)\ 


Im(z)’ 


pT\ 

Jo 


g{ia + Re'^) 


Re 


,i 0 


i?e'® + ia — zq 


dO 


< 


ll5lU=» [ 

Jo 


R 


0 |a + i? sin 0| |i?e*® + ia — zq 


■ de, 


which, by dominated convergence, converges to 0 as i? —>■ +oo when a is fixed. On the other 
hand, g{- + ia) belongs to while (•+!« — is in since ia — zq has a non-zero 

imaginary part. Therefore, the limit R —> -|-c» can be taken in the first integral on the right-hand 
side of (105), which leads to 


= & 


g{u} + ia) ^ \ / ■ \-i\ 


U! + la — Zq 


The conclusion now follows from the strong convergences of (• -I- ia — zq) ^ to (• — zq) ^ in (Ri^) 
and of g{- + ia) to g in as a —>■ 0. 


Proof of (iv). Let iphe a. real-valued function in From (12), we get 

J^g{uj + ig)(p{uj) duj = ^ J^{g,{- - uj - i??)"^)^-!^^! ^^(w) dw. 

Taking the real parts of both sides, we obtain 

/ Re (g{u! + ig)) <f{uj) dw 

Js. 


= 3^./ I 


— OJ 




+ ( Reg, ■ 


(106) 


(/?(w) dw. 




Consider now (j) G (^““(K^), with support contained in [—i?, i?] x R for some finite i? > 0. Then, 
Fubini’s theorem for distributions (see [38, Chapter IV, Theorem IV]) asserts that, for a given 
distribution T G o5^'(R) and tp G 


/ {T, (/){■, u})) pi^j) duj = (t, (j){-,ui)ip{uj)duj\ 
Jr ’ \ Jr / 
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When T G H ^(R), the above linear form can be extended to functions in Therefore, 

(106) can be rewritten as 


/ Re {g{uj+ ir]))ip{uj)duj = (img,^ [ — 

\ Jr (* 


-dw 

+ T / ff-i, ffi 


1 


iu) 7] / ff-i Ml 


(107) 


In view of the following strong convergences in i7^(M), 


-f 

27'" Jr 


1 




I 

Jr 


ip{uj)duj — 

V^o I 


the equality (107) leads, in the limit 77 —>■ 0+, to 




The first equality in the statement of item (iv) is finally obtained with the following lemma (recall 
that, according to Lemma 4, i7®(R) is stable by the Hilbert transform). The second equality 
follows by applying Sj to both sides and remembering that = —Id. 

Lemma 76. Let s ^ 0. For any T G i7“®(R) and Lp G i7®(R), 


ip)H^I — (2^;-^(7’))iJ-«,_f/o ■ 

Proof. Consider first the case when T,p G (R). Then, using Plancherel’s formula, the duality 
product can be rewritten using a L^-scalar product 

{S)T,p),^p^ = (Wr,4,)L^=2^ (j--i(.^r),J-V)^, =27r(-isgn(-).F-ir,J-V)i2 
= 27r(j'“iT,isgn(-)J^"V)i2 = - = -{T,^p)^py. 

The conclusion is obtained by a density argument. □ 


6.3 Proof of Proposition 12 

The proof presented in Section 6.2 can be followed mutatis mutandis upon introducing, for given 
elements f,g G H, the bounded causal function 


a/.9(^) = {f\Tc{T)\g) , 

and noting that ||a/,g||L=o ||Tc||L=(eCH))ll/llll5ll- 

The only additional technical point is the strong analyticity property, which is however easily 
obtained from the following bound: for z = w + ir; G U, 


dre(z) 

dz 


B{n) 


roo 

/ l’c(T)(ir) 

Jo 




Bin) 


^ \\T, 


c||L~(R^.Z3(H)) 


f 

Jo 


re dr = 


||'rc||L“(R^,73(«)) 


< +00. 


6.4 Proof of Proposition 13 

For z G U, we have 

4 c ( 2 ) = [ Ac{t)i 
Jr 


’ dr = —i 


c +00 

'/ ' 
Jo 


... — [tH Azt , 
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A simple computation shows that 

Ic(w + i? 7 ) = -i r ! dPf dr = / - - -- dPf = (w + ir; - i7)-\ 

Jo Js. Js. ^ + ’'■V ~ ^ 

The series of equalities can be made rigorous by testing them against functions f,g € H, and 
using Fubini’s theorem to justify the exchange in the order of integration. 

The strong convergence of Ac(- + i??) to in is ensured by Proposition 12. The 

Fourier transform can therefore be deduced from this limiting procedure. We consider the limit 
of Im Ac(- + iry), the real part of Ac{- + it]) being obtained from (14) and Definition 8 . 

Let f GH and ip G Then, using Fubini’s theorem, 


ImAc(- + ir?) 






II 

JR JR 


where the measure is defined by fij (b) 


(^_A)2+^2 <^MMf(dA)d(u = - 
= (/ l^b^l /) for any b G ^(R), and 


trjW (dA), 
(108) 



_n _ 

A)2 _|_ 


ip{uj) dw 



(/5(A + 7?C) d^. 


Note that 


|t^(A) -7r(/3(A)| 


^2 + 1 


I<P(A + -qC) - <p(A)| d^ < 


1 + ^^ 


dC, 


where the last bound is obtained by rewriting Lp[\ + q^) — <^(A) as the integral of its derivative 
and using a Cauchy-Schwarz inequality. This also shows that is uniformly bounded as ry —>■ 0+. 
Since the measure is finite, (108) leads by dominated convergence to 


ImAc(- +iq) ,ip 




->• —TT 


<P(A) /if (dA), 


which shows that ImAc = 


6.5 Proof of Lemma 14 


Let us first assume that ImTc ^ 0. The aim is to prove that ReTc ^ 0 on (—oo,a;o]. Consider to 
this end tp G oS^(K) with Supp((/ 5 ) C (—cxDjWq] and ip ^ 0. Then, for any uj ^ ujq and w' ^ 0, it 
holds (p{(jj — uj') = 0 , so that 


Vw ^ Wo, {Sjip){uj) = lim / 

e^ 0 + 


/R\[-e,E] 

Let f G TL. In view of (14) and Lemma 76, 


ii„ 

£—)-0 + 


r 


(109) 


/ 


ReTr 


= ((/ 


Imf; f) ,S)(p 






ImTt, 








The latter quantity is non-negative since ImTc ^ 0 and iji/? ^ 0 on Supp(ImTc) C [wo,-l-c») 
(by (109)). 

Let us now assume that ReTc 0 on (—(X),a;o]. The aim is to prove that ImTc ^ 0 on the 
support of this distribution, which is included in [(Wo,-|-oo). Consider therefore p G ^(R) with 
Supp((p) C [wq, -I-oo) and ip ^ 0. Note that 

Vw < Wo, (.^(p)(w) = lim [ ^ ^ dw' < 0, 

^_vn— / i.y 
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and, for any f € H, 


ImTr 




ReTc 


^0- 


This gives the desired conclusion. 


6.6 Proof of Lemma 21 

The fact that Bi G S(L^(K^),L^(K.^)) is a simple consequence of the following inequality: for 
ip G it holds, for almost all r G R^, 




R3 


Biir, r')(p(r')dr' 


^ l|-Sl(r, •)||i=o(R3) ||(/?||il(R3). 


This shows that Blip G L^(]R^) with 


\Biip\\L^R3) < \\Bi{r, •)||ioo(R 3 ) dr^ ||v 5 ||li(r 3 )- 


Now, for / G it is easy to see that Bif is an integral operator with kernel Bi{r, r')/(r'). 

In addition. 


\B 


i/lle2(L' 


/R3 JR3 


\Bi{r,r 


')/(r')|^drdr'< / [ ||Si(r, .)||io.(R 3 )|/(r')|^ drdr' 

Jr3 Jr3 


'R3 


■)llL2(R3)dr 1 ||/||i2(R3). 


This gives the claimed result. 


6.7 Proof of Theorem 24 

Fix 0 < 77 < a; and f,gG'H. We start from (24), which we rewrite as 
_ i r+°° - -- 

C+(z 2 + ia;) = — / A+[v + v'— ui'+ i{uj — rfj) Q B~ [v'— ui'— irf) Aui'. (HO) 

J — 00 

By Proposition 13 and Proposition 17, 

A+{z)=A*i{z-A2)-^Ai, JF{z)=B*i{z + B2)-^Bi. (Ill) 

The poles of z !->■ 4+(j/+z/'+i(w— 77 ) —z) are located on the half-line i(a;— 77 )-|-(— 00, i'+i''—a), while 
those of z I—>■ B~{v' — \r] — z) are located on the half-line —\rj+{b + v', -l-oo). For any closed contour 
not enclosing any point of those two half-lines, the integral of A+ {y+v'+\{uj—rj)—z^QB~ {v'—\rj—z) 
on this contour vanishes. Let us choose the contour plotted in Figure 6 and evaluate the 
contributions of the left-hand side of 



A+(y + v' + i(a; 


77 ) - z)gB (H - i 77 - z)f 


dz = 0, 


( 112 ) 


on the various segments. Recall that we choose v < a-\-b and v' G {—b, a — v)^ so that v + v' — a < 
0 < v' + b. Let us also emphasize that the operators appearing in the integrand do not have 
singularities. 
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Figure 6: Contour used to compute the integral (112), with 77 > 0 small compared to w. Note 
that the condition iz + v' — a<0<i^' + b ensures that the central vertical part of the contour 
does not intersect the poles of the functions in the integrand. 


Let us first consider the part of the integral corresponding to the right side of the contour. 
Using (111), we obtain that for all uj' G [0, L], 


( Tr-^ A+[v + v'— L + \{ijj — Tj — ijj')')gB — L — i{uj'+ rj)) f 
Jo L 

L 


duj' 


11 ^ 1 / 


0 


Tr^ 


AT 


_ ^ f 

V + v' — L + i(a; — ij — u') — ^ ^ v' — L — i(w' + 77 ) + B 2 ^ 


Ahj' 


i 


e2{n,nb) 

L 


\\BigL 




1 


1 

ly + ly' — L + i(a; — g — to') — A 2 

B{n^) 

v' - L - + g) + B 2 


dw' 


B{Hb 


^ C'll/llwllsllw 


dw' 


T , , , = C'll/llwbllwylog 

L Jq oj' + g L 


L + g 


where we have used A 2 — {ly + ly') + L ^ a — {v + v') + L ^ L. Similar estimates can be stated 
for the upper, lower and left parts of the contour. For the upper part for instance, for which the 
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integration is performed from z = L + iL to z = iL, we get 

rO 


Tr-H 


/: 

^CWfWnhWn f 

Jo 


^CWfWMn 


= C\\f\\n\\9\\n 


A+(i/ + v' — uj' + i(w — rj — L))gB~{y' — uj' — i(77 + L))f 


duj' 


1 


1 

y d- y' — w'd- i(a; — g — L) — A 2 

Bin.) 

y' — io'd- i(T d-g) — B 2 


dw' 


Bi-Hi 


1 


1 


L + r] J u}' + a — {v + v') 


duj' 


log + a - {v + v')) - log {a - {v -\- v')) 


L + r] 

where we recall that a — {ly + v') > 0. We then take the limit L —>■ +oo, so that the contributions 
to the integral which are not on the imaginary axis or on the real axis vanish. We deduce that 


f 


Tr^ 


+ y' — oj' + i(w — ri))gB~{y' — w' — ig)f 


dw' 




Tr-H 


' —OO 


A+[y + y' + i{uj - r] + u)'))gB-{y' + i(w' - g))f 


duj' 


Tr-H 


A+{y + y'-\-\{u}-\-Lo')')gB {y'+ iLo')f 


duj'. 


In view of (110), we finally obtain that 


C+iv + iuj) = - 




(* + 00 


A+ (h + h' + i(a; + w')) Q B (y' + iw') duj'. 


(113) 


We next note that our choices for y, v' ensure that the expressions on both sides are analytic for 
all w > 0, and can be extended analytically to all w G K. Therefore, the above equality also holds 
true for w ^ 0. 

In a similar fashion, we prove that, for all w G R, 

_ 1 f+oo _ _ 

C~{v+ ioj) =—— / A~{y+ y'+ U!')') Q B+{v'-\-\ijj')duj'. 

J — OO 

This equality is established as for (113) by considering C~{y — iw) for w > 0 and evaluating the 
various parts of the left-hand side of 


Trn 


% 


A~[y + y' — i(a; — g) — z)gB+{y' d-ig — z)f 


dz = 0. 


The poles of the integrand are on the half-lines —i{uj — g) + {y + y' + a, -foo) and ig+{—oo, —b + y'). 
The conditions y > —{a + b) and —a — y < y' < b ensure that —& -|- h' < 0 < h -|- h' -|- a, so that 
the integrand has no singularity on the imaginary axis. 

Finally, since A+(r)0i?+(—r) = A~ {t)qB~ {—t) = 0 for t ^ 0, we can concatenate C+(H-|-ia;) 
and C~{y + iw), and obtain 

~ 1 /■+°° ~ ~ 

C(^y itx^) = —~— / Aiy -\~ y i(tx^ -t- cu )^ O B{y it^ ) dec . 

6.8 Proof of Corollary 25 

The proof is based on the representation (25) with the choice w = 0. Consider y < a + b. It holds 

__ -r r+°° 

C+{y) = C+(y) = - 


27r 

1 

27r 


/ +CO _ __ ___ 

A^{iy jy' io;') O B~{v' + iuj') 

-OO 

y y' — Ao — iuj' 

0 


' —OO 
r*+oo ■ 


y d-y'- A 2 - iuj' . 

y y - A2Y + {uj'Y \ 


y' + B2- icu' 

^ {y'+ B 2 Y + {uj'Y \ 


dJ. 
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The odd terms in ui' cancel out by symmetry, so that 




^1 /,, I ../ /t \2 I ./'i2 


(Z2 + I/' - ^ 2 )^ + (w')^ 


© 


v' + B 2 

+ B 2 Y + {u:'Y ^ 


dui' 


At 


UJ 


{v + v' - ^ 2)2 + (w') 


/'i2 


0 




o; 


(!/' + B 2 Y + (wO 




(114) 


dw'. 


This shows that this operator is positive and self-adjoint in view of Lemmas 22 and 23. As a 
result, lmC'+(z/) = 0 for v < a + b. This proves the first assertion in (29). Also, we get from (114) 
that C'+(j 2 ) = Re (7+ ^ 0 for v < a + b. Together with Lemma 14, this shows the first assertion 
of (28). The results concerning lmC“ are proved in a similar way. 


6.9 Proof of Proposition 26 

The first assertion follows from the fact that the domain of Hm is iL^(K^^) and that TLatRc (4'^) = 
A^Re {d>%) and HNlm {d>%) = A^Im 

The density is bounded since it decreases exponentially fast away from the nuclei and is 
continuous [13]. 

In order to prove (36), we rely on (31) and (34) in order to write (recall thatis real valued) 

= n[ ([ g{r)d>%{r,ri,...,rN_i)dr] ( [ /(r')4'^(r', ri,..., rAr_i) dr'^ dri... drjv- 

jR3{iV-l) \J^3 J \Jk3 J 

= [ [ S(i')7w(i’:r')/(r')drdr'= (/jy^jg). 

Jr3 


To bound the kernel | 7 ^(r, r')p, we write 

l7w(i’.i’')r = [ ■■■ [ 4'w(r,r2,...rAr)4'^(r',r2,...rAr)dr2...drAr 

... [ |4'^(r,r2,.. .r7v)f dr 2 .. .drw) f/" ...[ |4'^(r', r 2 , • ■ • rw)f dr 2 . 

\Jr3 Jr3 / V7r3 7r3 

< PwWPw (!■'), 


where we used the Cauchy-Schwarz inequality. 

Let us finally recall why 2 defines a bounded integral operator. Note first that p% 2 (i'j I'O ^ 0- 
For f,g € Hi, the Cauchy-Schwarz inequality then leads to 


) dr dr' 

\ 1/2 


\{f\PN.2\9)\ = [ [ /(i')PAr,2(r.r')5(i’') 

^([ [ |/(r)|V^,2(r,r')drdr'y'7/ [ |5(r')|X,2(r, r') dr dr') 

\Jr3 7r3 / \Jr3 7r3 / 


\ 1/2 


(AT-l) 


^ (^-1) II ©I 
^ o \\Pn\ 


\ frp % 


1/2 


\9\^P% 


1/2 


l-Hillffllwi- 


This shows that p% 2 defines a bounded operator on Hi, with operator norm lower or equal to 

(^- l )|| p ^ L,=o / 2 .’ 


. .drw 
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6.10 Proof of Theorem 36 

Since 


A_/(ri,...,rAr-i) = ViV f /(r)^'^(r, ri,..., rAr_i) dr, 


and introducing 


N-l 


^N-1 = ^ Api, 


it is easily seen that (1—A^v-i)^- is an integral operator with kernel [(1 — Atv-i)^^^] (ri,..., rw-i; r). 
As G it follows that (1 — A 7 v-i)A_ G 62 {'Hi,'Hn-i)- Therefore, any operator of 

the form A*_BA-, where the operator B on 'Hn-i is such that (1 — AAr-i)~^'^^i?(l — G 

is trace-class. In particular, the operator 

a.Gh(T)|^^g_ = -A*_{Hm_, - E%)A_, 

is trace-class. 

Let us now compute more explicitly the action of this operator. Let 

•— 2^ ^ext- 


We use the definition (31) of a(/), and obtain 


^ /ii(ri)^ [o(/)l'w] (J*!,... ,r7v-i) = y/N /(rAr) hi(ri)^ . ,rAr)dr7v 

= Vn[ /(rAf) (iLo.w^w) (J*!, • ■ ■ ,rAr) drw (115) 

iR3 


[ (hi/)(rAr)l'?^(ri,...,rw) 


dr 


N, 


SO that 

(ilw-iA-/) (ri,..., rw-i) = E% (A_/) (ri,..., rw-i) - (A_/ii/) (ri,..., r^-i) 

Moreover, it is easily seen that, for any G Hn-i, 

{A*_^n-i) {r) = VN / 4'^(r,ri,... ,rAr-i)$Ar-i(ri, ■ • • ,rw-i)dri .. .drw-i. 

1r3(N-1) 

{A*_ {Hn -1 - E%) A_/) (r) = - {j%hif) (r) - f Kn^t, r')/(r') dr', 

iR3 

KNirr')=N^ [ ^^(r, n, ■ ■ ■ , rAr-i)l'^(r', ri,..., r^-i) 

■“ yR3(iv-i) 


Therefore 


with 


r' - r,; 


■dri .. .drAT-i. 


Since we already know that the integral operator on TLi, with kernel KN{r,r') is trace-class, 
and that 'i’% is continuous and decays exponentially fast (see e.g. [13]), we have [40, Theorem A.2] 


Tr«, (Kn) = [ Kw(r, r) dr = N^N - 1) [ 

Jr3 Jr3(N-1) 


l'^(r, ri,...,rAr-i)| 


F-ril 


■ dri... dr^v-i 


r f 

W |r-r'| 
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where we recall that 2 i® the two-body density matrix defined in (35). 
Finally, 


Tr^, (5rGh(T)|^^„_) = -Tr^, - E%)A_) 


= Tr^i {l%h^) + 2 f f ^ dr dr' 

Js.3 7 r 3 |r — r I 


which gives the claimed result in view of the following representation of the ground state energy: 

Pn,2(^a') 


E% = = Tr^, {h^^%) + f f 

iR3 Jr3 


dr dr'. 


(116) 


6.11 Proof of Lemma 42 

Proposition 26 shows that p% G LP(M.^) for 1 < p ^ -l-oo. This implies that '^n 

belongs to Hn for u G C since, by Holder’s inequality, and the inequality ^ 


NJ2i=i ’vi'Tir, it holds that 
2 


[ \'i’%(ri,...,rN)fdri...drN^Nf p% ^ N\\v\\14p%\\l3/2 

Jm3N J J^3 

^ NCc4v\\1,\\p%\\^3/2, 

where we have used the embedding C' ^ L®(]R^). Moreover, p'jf G L®/®(IR.^) ^ C, so that 
\{v, p%)c',c\ ^ Ikllc'lIPwllc- We therefore deduce that H is a bounded operator from C' to Hn, 
whose norm satisfies 


\B 


\b(C',Hn) ^ y^7VC'c'||p^||i3/2 -I- ||p?/||c- 


We finally have, for v G C', 


K 


N 


^(g) 


'^n) =[ v{r)p°j^{r)dr = {v,p%)c',c 

In. 


from which we deduce that {^^4Bv)-h^ = 0. Since v was arbitrary, we conclude that 5*4^ = 0. 


6.12 Proof of Theorem 46 

Consider f,g £ C'()°(R^,C) (that is C°° with compact supports). Using the fact that {Hpf — 
i?^)|4'^) = 0, we obtain 


{f,v-^B*{HM - E%)Bg)c',c = (/ \B*{Hn - E%)B\g)c' = {Bf \Hn - E%\Bg)^^ 




(Hn - E\ 




We next observe that 


{Hn-E%) ^5^p(r,)4'^(ri,...,rw)^ = 



N \ 

5]] A5(ri) 4'^(ri,.. 


2=1 


N 

^ n ) - ^5(u) ■ Vr,4'^(ri,..., rAf), 

2=1 
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so that, using the fact that is real-valued (see Proposition 26), 


{f\B*{HN-E%)B\g) =^2 f /(rj) l«'w(i'i.---i‘Ar)pdri...driv 

ij=i ^ 2 

^ r _ 

- y" / /(rj)^'^(ri,...,rAr)Vg(r,) • Vr,^'^(ri,...,rAr)dri...dr7v 

2 r _ 

= 1 ^'^ V/(rj)-V5(r,)|4'^(ri,...rAr)pdri...drAr 

+ / /(rj)V5(r,)-Vr, (|«'w(ri,---,riv)|^) dri...drAr 

N r _ 

- y" / /(rj)^'^(ri,...,rAr)Vg(r,) • Vr,^'^(ri,...,rAr)dri...drAr 


V/(r) • V 5 (r)p 5 ,,(r)dr. 


We deduce that 2v^ ^B*{Hpf — E^)B = 2v^ ^B*{E[‘2 — E^)B = —div (p^V-) as operators on 
the core C'“(R^,C). We next observe that div (p^V •) can be extended as a bounded operator 
from C to C. Indeed, for f,g G C), 


( 5 ,div (pwV/))c,_c ^ Ik^lLoo I|V/||l 2 ||V 5 ||l 2 = 47 r||p^||Loo||/|lc/|| 5 ||c', 
which shows that 

||div (pwV-)||b(c-,c) < 47 r||p^||^^. 

Therefore, 2v~^B*{E[2 — E^)B = —div (p^V-) as bounded operators from C to C. 

For the second part of the proof, we first note that it is sufficient to check the convergence in 
the case when f = g G C'“(R^,R). It holds: 


f,v:^B* 


Hi - E° 


N 


- {f,v-^B*{Hi - E%)Bf)c',c 

/ C',C 

{h 2 - E%r 


iH2-Ei)^+co^ 
= -(j,v-^B* ' 


iHi-E%r+uj\ 


Bf 


C',C 


r+co 

Jo 


A3 


0 


ttU tpQ 

pH^ 


'Hn 


where we used the spectral decomposition of — E^ in the last equality. The integrand of 
the last integral converges pointwise to 0 when |w| —>■ + 00 . It is also non-negative and uniformly 
bounded by Al[o^+oo)(A), which is integrable since 


r +00 

/ Ad 

Jo 


P 


’(Bf) 


T-Ln 




< 00 . 


The weak convergence therefore follows from the dominated convergence theorem. 
To prove the strong convergence, we use the following rewriting for g G C: 


v-^B 




{H2-E%r+co\ 


Bg = v-^B*A^Mg, 


where 


A.U, — 


(Hj - E%f 
{Hi-E%Y+u:^ 
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strongly converge to 0 on 'Hat, and 

^ / 1 \ 

Mg := - E%)Bg = ^ --Ag(r,)vl/^ - V5(r,) • . 

When g G C' is such that Ag G LF‘, it holds that Mg G TLn (by a proof similar to the one in 
Section 6.11), which allows us to conclude. 


6.13 Proof of Proposition 48 

We first prove ( 68 ) and (69). We start from an expression similar to the one provided by Lemma 38: 

Go{z) := 4o,+ (z — (ilo,w+i ~ ^o,n)) ^o,+ + ^o,- ~ ^o,n) 

Then, we notice that, for f G Hi, it holds 

N+l N+1 

hi{ri) = Y ^i(i’*) («^(/)l^w)) = aHhif) \<^%)+a'>(f) 

Z—1 Z —1 

or, equivalently, 

Ho,N+iAl+if) = AI+ (hif) + E%Al+ify, 

so that 

{z - iHo,N+i - El^)) Al^ = Al^iz - hi). (118) 

Hence, the particle part of ( 68 ) is a consequence of the equality = 1 -^ — 7 ^ ^ (similarly 

to (39)). To handle the hole part, we use computations similar to (115), and find 

lto,iV-i4o,-(/) = EQ pfAo-{f) — Aq- (hif) . 

We deduce that 

(z + Ho,n-i - < 7 v)^o.- = Ao,-(z - hi), (119) 

and we conclude using the fact that 4g _4o,_ = Jqjy- Combining (118) and (119) leads to 
Go(z)(z — hi) = Aq _4o,_ + 4o,+4q _|_ = Upon replacing z by z and passing to adjoints, it 
also follows (z — hi)Go{z) = This shows that Go{z) = {z — hi)~^. 

To prove the first assertion of Proposition 48, we notice that the operator-valued functions 
T !->• —i0(T) (l-Hj — 7o A/) and t >->■ G'o_p(r) have the same Laplace transforms (see Propo¬ 

sition 13). We conclude that the two operators coincide since the Laplace transform is one-to-one. 
The proof for Go,h is similar. 


[hi{n) 




(117) 


6.14 Proof of Proposition 50 

Let z G C \ R. Let us first prove that G(z) is a one-to-one operator. Let f G Hi he such that 
G(z)/ = 0. From 


Im 


{{f\G{z) 


/) =Im 


“HAr+i 


(^(^Alf\{z-{HN+i-E%)) 'I 4 ;/) 

+ (4_/ (z - {E% - Hn-i))~" A_f) ) 

= -Im(z) [ ((Re(z) + U^-A)2+Im(z)2)-^d 
Jr 

-Im(z) [ ((Re(z)-U^ + A)2+Im(z)2)-^d 

Jr 


( 120 ) 


P. 


Hn+1 


P- 


Hn-1 


(Kf) 

{A-f) 


'Hn+1 

2 


'Hn-1 


( 121 ) 
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we deduce that both terms in the right-hand side must vanish. In particular, since Im (z) ^ 0, it 
must hold A’^f = A-f = 0. In view of the identity A*_A- + A+A*\_ = this implies / = 0. 
Hence, G{z) is one-to-one. 

As a consequence, G{z) is an invertible operator from 'Hi to its image D{z). Since ^G'(z)^ = 
G{z) is also one-to-one, D{z) is dense in Hi. 

Let us finally prove that D{z) C We use to this end the equality (51). Let us consider 

the hrst term in this equality. A simple computation shows that, for any G Hn+i, 

(A+$7v+i) (r) = VN+1 / $Ar+i(r,ri,..., rAr)^'^(ri,... ,r7v)dri. ..dr^, 

Jr3N 


so that A+ is a bounded operator from H'^{'E.^^)r\HN to Since A’^ is a bounded operator 

from Hi to Hn+i and {z — H^+i + is a bounded operator from Hn+i to C^Hn, 

we deduce that A+{z — Hff+i + E^)~^A'^ is a bounded operator from Hi to Similarly, 

for any ^n-i £ Hn-i, 


{A*_<i>N-i) {r) = VN [ 

Jr3(N-3) 


$Ar_i(r 2 ,..., r7v)^'^(r, r 2 ,..., tat) dr 2 ... drAr, 


so that A't_ is a bounded operator from Hn-i to This allows us to prove that Al(z -|- 

Hm-i — is a bounded operator from Hi to Finally, G(z) is a bounded operator 

from Hi to which proves that D{z) C 


6.15 Proof of Lemma 56 

Let us first prove that defines a bounded causal operator. The proof is similar for 

rely on the following result. 

Lemma 77. For all h G the operator Jq jifh is a Hilbert-Schmidt operator on Hi, and 

there exists K G such that 


yhGL%R^), ||7oVIIs,(h,) 

Proof of Lemma 77. Since 7 ° is a projector, for h G L®(R^), the operator 

%o.w7o,JV^ = Ml - A)-^/Ml - A)^/Mo,7v( 1 - A)i/Ml - 

is the composition of (1 — A)^/^7g ^(1 — A)^/^ G 6i{Hi) with the bounded operators (1 — A)“^/^/i 
and h{l — A)“^G_ jn fact (1 — A)“^/^h G &e(Hi) with 


( 1 _ A)-i/ 2 /i 


eeiUi) 


^ K\\h\\L6 


by the Kato-Seiler-Simon inequality [39, 40]. Therefore, Jq jyh G 62(7^1) with 




esCHi) 


< A:||h||i 6 , 


□ 

, C), it holds 
( 122 ) 


which concludes the proof. 

We now proceed to the proof of Lemma 56. We first note that, for f,gGG[ 

= {Vc^^^f\9)c = {vV^l,g)c’,c- 
In particular, for t G and for f,gGHi,we get 

(/lM°ym(Mff)-Hi = (/ = -i 0 (M (vy^J,Go,p{T) 0 Go.h(-T)wy^g^^^ 

= -i0(r)Tr7^, Go,p(r) (vl^^g^ Go,h(-T) 
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Let us prove that G'o,p(t) Go,h(—T) is indeed a trace-class operator. Replacing 

Go,p and Go,h by their expressions found in Proposition 48, and owing to the fact that 7 q ^ is a 
projector commuting with hi, we obtain 


|(/l^symMl5)| = Tr«, -7o.w)e 7o,Ne''"'''7o,N7o,N f)) 


0 ir/ii 


^ II ~ 7o,n)^ 




IS(Wi) 




B{Hi) 


According to Lemma 41, G C ^ L®(IR^). Therefore, 7o G by Lemma 77, 

hence is bounded, with 


7o.7V 


,1/2 




Bi-Hi) 




7o,n 




S 2 CH 1 ) 


^G 


yv < cii/ii«, 


LO 


Similarly, 


(vl^^g) 7 S,n 


of r such that 


BCHi) 


^C/||5||hi- Altogether, we found a constant G G M"*" independent 


Vf,gG Ui, |(/|Py+(r)|g)| ^ G\\f\\HM\n,, 
which proves that is a bounded causal operator on 'Hi. 


Let us now prove that 4>k is a bounded operator from Hi to C. Recall that 4>k is real-valued 
and (t)k G H2(k3) ^ pi ioi l^ks^N. For / e Hi, we obtain 

ll'/'fc/llc^C^II'/'fc/llLa/a 110^11^3 I|/||w„ 

where C is a constant independent of /. The proof that (pk is also a bounded operator from C 
to Hi is similar, noticing that C ^ L®(R^). We now use (80), and find that, for f,g Ghii and 
for r G K+, 


(/l-PsymMIff) = -7o.w)e |^fc)e‘^'''“(<?ife| (i^yV)) 

k^l 
N 

k^l 

N 

= -''^'^{{yl''^f)4>k (l-Hi -7S.Ar)e"‘^^'''"'''“yi'Hi -7S.Ar) {^1'^g) (pk) , (123) 


fc=l 


which gives (81). 

We finally prove that H,y’“(r) = H!1+(—r). Performing similar calculation as for we find 

that, for f,gG'Hi and for r < 0, 

N 

{f\Psy^{r)\g) = -iJ2 ((I'y^^) </'fc (l-Hi - 7S,Ar)e'^^^'“®'“^ ^t^k) 

k^l 

= (5l^’sym(-'r)|7)- 


For a bounded operator A G H(Hi) and for f,g G Hi, it holds {f\Ag) = {Ag\f) = {g\A*f), so 
that, since the functions pk are real-valued for 1 ^ k ^ N, 


(/|H°-(r)|5) = (f\{P^U-r)y g) ■ 


(124) 


57 






































Since hi is real-valued, in the sense that hif is real-valued whenever / is real-valued, we easily 
get that _ 


V/ G Hi, {In, - 7oV)e-'"'‘V = (Iwi 




so that. 



Together with (124), this proves P°’ (r) = P°’+(—r)*. 


6.16 Proof of Proposition 59 

The expression for in (84) comes from the expression for in (82). Since for k ^ N, \i 
holds £k ^ £n, we obtain 

yi^k^N, {t-Hi - - £k) > £n+i - £n > 0- (125) 

From (84), we deduce that PgO,jj(ia;) is a negative bounded operator for all w G The self¬ 
adjointness comes from Remark 57. 


The bound (85) is proved similarly as in Lemma 29. Let us now prove (86). From (125), it 
holds that, for all 1 ^ fc ^ A^, 


0 ^ (lwi~7o.Af); 


hi - £k 


^ sup 


E 


+ {hi - Skr E^elZ-e. W + 

In particular, there exists a constant C G M"*" such that 


1 

2w 


if w ^ ew-i-i — £n 


£n+i — £n 


up -f (eAT+i — EAf)^ 


otherwise. 


Vl^A^iV, VwGR,., 0^(1«,-7o%) 2 a 

-I- (hi - EkY Va;2 -|-1 


Using the fact that for 1 < h < A^, (j)k is real-valued, with J2k=i — Po Af> obtain 


Vw G 


0 < -Ps%(iw) < 


2C 


N 




-bi 




= 


2C 


k=l 


\fijp 




which proves (86). The fact that is indeed a bounded self-adjoint operator on Hi 

comes from the fact that vY is a bounded operator from C to Hi and from Hi to C, and that the 
operator of multiplication by (j)k is a bounded operator from C to Hi and from Hi to C. Together 
with the fact that the function oj i—>■ (w^ -I- 1)“^/^ is in LP{R,^) for all p > 1, this implies that 
Isym(i’) € L^{^uj,S{Hi)) for all p > 1. The analyticity of this map is straightforward. 


6.17 Proof of Theorem 60 

Let us prove the equality 2 4>k{^'Hi ~ 7o w)(^i ~ £k)4>k = —div (pj) , as operators from 

C to C. We first note that, since (pk G P"‘(R^) for 1 ^ h ^ A^, it holds 4>k(t>i ^ for I ^ k,l ^ N. 
In particular, 


N N N 

E= E E \pk(pl){£l — £k){(pl(pk\ = 0 , 

k—1 k—11—1 


SO that 


N 

2 E Pk(t-m - lo,N){hi - £k)pk 


N 

2 ^ ^ (j^kihl ^k)4^k- 


(126) 
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Consider now f,g G C). In view of the equality 


{hi - £k) {(j)kg) = (l)k ( “2^5 ) - V^fc • Vg, 


it follows that 

N 


2 / 


^ ^ 0/c {hi ^k)4‘k 


k=l 


9) = / 


Hi 


N 


{-^g)) - 2 (/ 


k^l 


ni 


N 




fc=i 


Hi 


{f\p%{-Ag))ni - (/|Vp° ^ • Vg)ni = [ p^n 

Jr3 


V/ • V5, 


where we used an integration by part to obtain the last equality. Together with (126), we ob¬ 
tain that the operators 2 J2k=i 4’k{^Hi — 7o Ar)(^i ~ Sk)4’k and —div (pj] ^V-) are equal on the 
core C“(K^, C). The end of the proof is similar to the one of Theorem 46. 


6.18 Proof of Lemma 62 


It is sufficient to prove that, for any v' < £n+i, v > £n — a;,a;' G Ki^, and f,g G Hi, it holds 

Tr-Hi Go,h{i' Hi'' + i(w -h w')) GoA’^' + 

= (GoA^ + 1^' + i(w -h w')) (vlGg^ GqA^' + 

We will only consider the case ly = 0, v' = po, to = uj' = 0 ioi simplicity, the other cases 
being similar. We rely on the fact that if A, B G B{Hi) are such that AB and BA are trace- 
class operators, then Tr(Ai?) = Ti'{BA) [40]. In our case, we consider f,g G C Hi, so that 
fi := and gi := vl^'^g are in C fl T°°, and we set A= fi and B = Go,h(Mo)<7iGo,p(/io). The 
operators A and B are bounded operators on Hi. Moreover, from Proposition (48), we get 

BA = GAMgiG^MTi = 7oV ~ A 

’ \po -hi J \ po-hi 



It holds 7 Qjy G 6i(Hi). Also, from the definition of po, it is easy to see that there exists 
0 < c < G < oo such that 

c(l - A) ^ \po - hi\ < G(1 - A). 

In particular, the operator {po — hi)~^gi = [(po — ~ A)] [(1 — A)“^(jii] is the composition 

of two bounded operators. Actually, the operator (1 — A)“^pi is in the Schatten class &q{Hi), 
thanks to the Kato-Seiler-Simon inequality [39, 40], and it holds 

11(1 - A)-V||b(„,) < ||(1 - A)-V||s,(«,) ^ G\\gi\A ^ C'Mni, (127) 


where C G K'*' is a constant independent of g. Similarly, (1 — 'Jo,n){po — hi) ^/i is a bounded 
operator, satisfying estimates similar to (127). Altogether, we deduce that, for all f,gG G)]“(R^), 

GoApo) (vlGg'j GoApo) e 6i(^i), 


with 



epHi) 


^ C'IIsIIkiII/IIhi, 


(128) 


where G G K+ is a constant independent of / and g. The proof that AB = Go,p(/ro)5iGo,h(/^o)/i 
is trace-class is similar. As a result, we deduce that for any f,gG G“, 


Tr-Hi ((^yV) Go.h(Aio) [AAg^ GqApo)) = Tr«i (Go,h(Ato) {vAg'^ GqApo) (^c^V)) • 


The proof for the general case f,gGHi is deduced by density from the estimate (128). 
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6.19 Proof of Proposition 73 

Let us first prove that s is a bounded linear operator. Let G'®-PP(/ro + i-) S S('Hi)). For 

f,g€'Hi and w S 


gapp 


/ + 00 , _ . _ _ 

Tr-Hi [G^PP{tJ-o + i(w + uj'))gWP{iuj')f 


dw'. 


Let us first treat the exchange part K^, and prove that it is a Hilbert-Schmidt operator. From 
the definition (99), is an integral operator, and, from Proposition 26, its kernel satisfies 



\K, 


I /m 2 j A ! ^ f f ^o,w(’')Po,Ar(''0 j j , 

,(r,r)| dr dr < / / - , -dr dr 

|r — r I 


< oo. 


where the last inequality comes from the fact that Po * | • | ^ S (since Po ^ S L^(]R^) fl 

L°“(R3)^ while I • € L\R^) + L°“(R3)) and pj]_^ G L\R^). We conclude that is a Hilbert- 

Schmidt operator, hence is bounded. 


For the correlation part, we use the following lemma. 

Lemma 78. For all f,g € %! and all u G K(.j, the operator gW^{iu))f is trace-class, and 

C 


3CGR+, yf,gGni, 


gW°{iuj)f 


€ 


^ f 2 I TM/o J -Hi g Hi- 

SiCHi) + 

Proof of Lemma 78. We first prove (129) for f = g G Hi. Let if G C')?°(]R^, C), we have 


if 


fW0{iu;)f 


if) = (if 




if 


and we infer from (95) that 

V/ G Hi, Vw G R,., 




C 


PH,) ^ (w2 + 1)1/2 


(129) 


Since 


|/^^cPo.w'*^c/| 


Si(Wi 


Po^TV^c/ 


62(^1) 7r3 Jr3 


/ / 

./R3 JK- 


Po.AfW i/(i-')r 


r — r 


J\2 


drdr' < C||/||^^, 


where we used again the fact that Po at * | • | ^ is bounded, we obtain that (129) holds true for 

f = g. For f ^ g, we deduce from the fact that xOyjjj(ia;) is a bounded self-adjoint negative 
operator, that 


gWff{iuj)f 


epHi) 




eiCHi) 


ff^^yV-Xsvm(iw)V-Xs°vm(iw)^'yV 


61CH1) 




< 




1/2,/_^0 


Xsym(i^) 



^J-x^y^yyf 

62 (Wi) 



c 


1/2 


(w2 + 1)1/2' 


61CH1) 

•hiIIpII-hi- 




e2(Hi) 
1/2 


61CH1) 


□ 
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We now proceed to the proof of Proposition 73. From Lemma 78, we get, for f,g € Hi, 

+ 00 _ 


Gapp 


(/r + iw) 


1 

2tt 


' —OO 
p + oo 


TyHi (G'^PP(7‘o + i(w + w'))gfF°(iw')/j dw' 




(* + 00 


<iC’ 


1 - - 

-J G-PP{po + K^ + u;')) 

G^ino + K^ + Lo')) 
G^^{go + p) 


BCHi 


gW^iiLo')/ 

C 


L2(R„,BCHi)) 


B(Hi) (a;'2 + 1)1/2 
'Hi ll^/ll'Hi , 


doj' 

eiCHi) 

Hibll-Hidw' 


where we used Cauchy-Schwarz inequality for the last inequality and the fact that uj i—>■ + 

l)-i/2 g Here, G' does not depend on w G nor on f,g G Hi. Altogether, we deduce 

that s is a bounded linear operator from L'^{R,^,B{Hi)) to L°°{M.,^,B{Hi)). 

We now prove the claimed properties of flA- Consider M > 0. By definition of go, the real 
number d := max(eAr_|_i —/tq, /tq —Eat) is positive, and \go — hi \ ^ d. Let us choose 0 < Am < d/M. 

C M, 


For 0 < A ^ Am and (^o + h) G L°°{M.,^,B{Hi)) such that E'^pp(^ o + h) 

it holds 




Vuj G Kij, /To + iw — hi — AE'^PP(^o + iw) = [mo + iw — hi] ^1 — A [go + iw — hi] ^ E®-PP(/ro + iw 
Since 


X[go + iLO — hi] ^ E*‘PP(^o + iw) 


A 

s; -1 

sCHi) d 


E^PP(Aio +iw) 


^ < 1 , 

B{Hi) d d 


the operator 1 — A [^o + iw — hi] ^ S®-PP(/io + iw) is invertible, with 


(^1 - A [/Lto + iw - hi] ^ E®-PP(/ro + iw)^ 


-1 


€ 


B(Hi) ^ XmM 

Since go + iu) — hi is an invertible operator with ||(/ro + iw — ^ + d‘^y^'^, we obtain 


that go + iu! — hi — AE®'PP(/ro + iw) is invertible, with 

1 -1 


go + iui — hi — AE'^PP(/io + iw) 

We deduce from this inequality that 
OA (/^o + h) 




BiUi) XmM 


d iq , . , ,_i|| , Km 

[go + iuj hi) llecH^) < ^^2 1)1/2- 




+ 


0 A (e-pp) (go+i-) 


L2(R„,B(-Hi)) 


^ Cm, 


where the constant Cm G R'*" does not depend on A G [0, Am]- This gives the claimed result. 
Finally, (102) is a direct consequence of the resolvent formula. 

6.20 Proof of Theorem 74 

Let us denote for simplicity 


and fix M > ||s|| C°(g + i-) 

L2(R,,,S(-Hi)) 

choice of M > 0, and introduce 

M 


S|| — l|s||B(L2(R^,eCHi)),L“(R„,B('Hi)))! 

. Let Am and Cm be chosen as in Proposition 73 for this 


C°{g + i-) 


L2(R,..B(Wi)) 


> 0 . 
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For this choice of r, it holds that, for any G'^pp g S ^Gq, , 


G^pp 


Therefore, from Proposition 73, flA o s 


(Mo + i-) 
G^P 


M. 


is well-defined for all A G [ 0 , Am]- 


Let us prove that there exists A* > 0 sufficiently small such that for any 0 ^ A < A*, gA os 
maps f 8 ^G°,rj into itself. For G'^pp G S ^Go,r^, it holds 

-G])) (Aio + i-)| 

fflAos G^-PP -oaos Gc 


flA os 


L2(R„,Z3(Wi)) 


) {fJ-0 + i-) 




+ 


,0AOS 


Go 


— G° ) {fiQ + i-) 


L2(R^,eCHl)) 

To control the first term (130), we use the following result. 

Lemma 79. The map 0 a os is (AG|^ ||s||) -izpsc/izte on 55 ^Go,r^ . 

Proof of Lemma 79. Let G^^^, G^^^ G 55 ^Gq, r'j . From (102), we obtain 


(130) 

(131) 


0AOS G^PP -0AOS 

From Proposition 73, 

0A os 

Moreover, 


Gapp 


= A(0Aos[Gf5]) (, 


^app 


— 5 


^app 


9\0 5 


^app 


(132) 


^app 


(Mo + i-) 




(s[Gr]-s[^])(/co + i-) 


L°°(R^,l3(ni 

Plugging these estimates into (132), we obtain 


^ llsll 


< Cm for 1 ^ j < 2 . 


(Gapp-c-app) 




,0AOS 


G^pp 


- 0A os 


G^pp 


(Mo+i-) 

(^-^) (po+i-) 

which proves that gA o s is (AG|^||s||)-Lipschitz on 55 ^Gq, rj . 


^ AGmI|s|| 


L2(R^,i3(Wi)) 
L^(R^,BCHi)) ’ 


□ 


Let us now control (131). By noting that 0 a=oos Gq = Go, we get from the resolvent formula 


that 


gA os 


Gn 


Go — (0A ~ 0o) OS ^Go^ — A ^gA os Go ^ ^s Gq J Go- 


Using estimates similar to the ones used in the proof of Lemma 79, we deduce that 


.0AOS 


Gn 


~ Go) (po + i-) 


^AG^ 


Gn 


(to + i-) 


..-B(Hi)) 


L^(R^,BCHiy) 

From Lemma 79 and (133), we arrive at the conclusion that for all 0 ^ A ^ A*, where 
A:,: - 


(133) 


Gm I lls|k + 


Gc 


(To + i-) 


..-B(Wi)) 
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it holds 0 A os (Go,?-)) C » (Ccrj. 

Finally, without loss of generality, we can assume that A^Cl^HsH < 1, so that, from Lemma 79, 
we get that for all 0 ^ A ^ A*, the map 0 a os is a contraction. The end of the proof follows from 
Picard’s fixed point theorem. 
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